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ABSTRACT 


A detailed analytical and experimental study of the sensitivity of 
Wave Digital Filters has been conducted. The results indicate that the 
wave digital filter tends to achieve the same low sensitivity character- 
istic as the analogue circuit from which it was derived. Other results 
indicate relatively higher sensitivity to terminating resistance values 
compared with internal element values, lower sensitivity for algorithms 
derived from simple rather than multiple elements, and higher sensitivity 
at the critical frequencies. Finally the rms error due to the quantiza- 
tion in the number of bits in the multiplier coefficients has been 


measured at approximately 3 db per bit for the many examples tested. 
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I. INTRODUCTION AND SCOPE 


A. INTRODUCTION 

Signal processing is necessary in diverse areas of science and engi- 
neering, such as communication, social sciences, biomedical, control, 
radar, acoustics, telemetry, and intelligence and information gathering, 
etc. In general this process can be done with analogue (continuous 
time), or digital (discrete time, discrete amplitude) signals. 

With the advent of LSI and VLSI, microprocessors, the need for effi- 
cient digital signal processing algorithms becomes more and more impor- 
tant. Much of the current literature is devoted to design of the linear 
algorithms, under the title of digital filters. Important factors in 
designing the digital algorithms are, time of the calculation, imple- 
mentation, accuracy (error), etc. 

Although there is evidence that direct digital filter design is 
possible [1], nearly all of the digital filter design algorithms use 
the analogue to digital transformation techniques. It is interesting 
to note that the rapid development of the digital signal processing is 
partly because of the existence of the well established theory on the 
analogue techniques and partly because of the abundance of the general 
purpose digital computers. 

The digital computer and specially microprocessors, being physical 
objects from space and economical point of view, can only accommodate 
for finite precision in the size of the digital filter multiplier co- 
efficients. Thus the need for digital filter algorithms with low sen- 


Sitivity to digital filter multiplier coefficients arises. 


17 








Fettweis [2] in order to design a low sensitivity digital algorithm 
has proposed an alternative digital filter design method, namely wave 
digital filter, in which analogue LC circuit is transformed into dis- 
crete algorithm using wave or scattering matrix parameters. This ap- 
proach in design is different from the conventional design techniques, 
because we have a new set of variables which are referred to as "incident 
and reflected wave parameters”. 

Wave digital filter design is believed to be difficult to understand, 
and is therefore generally avoided. But in reality this is not the case. 
In fact, to design a wave digital filter one need not go into the details 
of the algorithm development. One merely needs to know some basic facts 
and then can use the final wave digital filter equations and tables in 
ιο desio the required filter. 

It is a well known fact that analogue LC networks have very low sen- 
sitivity to variation in LC component values. Upon this fact Fettweis 
and others have argued that, since the wave digital multiplier coeffi- 
cients are derived from the LC parameters of the parent circuits, they 
should also have the same favorable low sensitivity characteristics to 
multiplier coefficients. Furthermore it is also known that the digital 
algorithms with low sensitivity to multiplier variations also exhibit 
minimum round-off noise due to quantization after multiplication of these 
multipliers with signals. As a result it is conjectured that the wave 
digital filter will have minimum round-off noise properties when compared 
with other digital filter algorithms. The purpose of this thesis 1s 


to analyze and check this conjecture. 
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B. SCOPE AND ORGANIZATION 

In order to develop and study the wave digital filter, the research 
in this thesis is divided into seven chapters. In Chapter II a brief 
discussion of the general digital filter theory is given. The presenta- 
tion contains only selected and necessary background required for wave 
digital filter theory development in this thesis later on. Also in- 
cluded in this chapter are A B C D matrix theory, and the concept of 
the delay free feedback which plays an essential role in the wave digi- 
tal filter theory development. Necessary sensitivity theory background used 
in the sensitivity analysis of digital filter is discussed briefly in 
Chapter III. The development of wave digital filter theory is done in 
Chapter IV. This development is straightforward and general in a sense 
that only one algorithm is developed for both series and shunt element. 
The effect of sampling interval is also introduced for the first time 
into the wave digital filter algorithms. Four useful tables of wave 
digital iterative algorithms for simple L and C elements in both series 
and shunt configuration are also given. The delay free path which plays an 
essential role in wave digital filter theory is emphasized throughout 
this chapter. The studies of the sensitivity of the wave digital filter 
to quantization in the number of bits of the multiplier coefficients is 
done in Chapter V. To do these studies in a fairly general sense, three 
different wave digital filter algorithms developed in Chapter IV are 
used with two different conventional digital filters for the given 
filter. A total of nine different filter types with different terminat- 
ing source resistances were examined and compared with each other to 
arrive at a general conclusion. Chapter VI presents a specialized in- 


depth study of the sensitivity of the simple wave digital filter 
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algorithm. In this chapter sensitivity distribution along the filter 
Structure was examined in order to understand the sensitivity behavior 
of the subsections of the wave digital filter on the overall sensitivity 
of the wave digital filter. Chapter VII summarizes the new results and 
proposes future research ideas related to the wave digital filter theory. 
There are three appendices. Appendix 1 includes an example to show 
the nearly exact equivalence between the wave digital filter and the 
conventional digital filter, both in the time domain and frequency domain 
when infinite precision arithmetic is used. The rest of the appendices 
include ten computer programs. The computer programs are used to derive 
the results of the main text. It is important to note that in order to 
facilitate the better understanding of the computer programs, explana- 
tory remarks are made in the comment cards. These computer programs are 
in FORTRAN IV and can be used on any standard general purpose digital 
computer capable of dealing with FORTRAN IV scientific computer language. 
Finally it is important to note that in this thesis for easy access, 
the references are given at the end of each chapter, rather than 


collectively at the end of the thesis. 
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II. GENERAL BACKGROUND 


A. INTRODUCTION 

The main intent of this chapter is to briefly review general digital 
filter theory in order to establish the background necessary for 
the main subject of this thesis, i.e. wave digital filters. The design 
of the conventional digital filter is a well established subject. Thus 
it will be dealt very briefly. Also discussed in this chapter is 
the A B C D parameter matrix theory, which will be used later on followed 
by the concept of the delay-free feedback path or delay-free loop in the 
digital two port network which are used in deriving the causal wave 
digital filter algorithms in Chapter IV. z domain transform theory is 


assumed as a background and is not discussed. 


BEEGENERAT REVIEW OF DIGITAL FILTER DESIGN 

The design of electronic filters in the analogue domain is a well 
established subject, which not only includes very sophisticated tech- 
niques, but also has some very well established supporting comuter pro- 
grams as well. Much of the development of the digital filter has been 
directed towards the transfer o£ these results from the analogue domain 


to the digital domain. 


In the most general sense a digital filter is a linear, shift in- 
variant discrete time system which is realized using finite precision 
arithmetic. The design of the digital filters involves three basic 
steps: 


(1) the specification of the desired properties of the system 


ID 
(~) 








(2) the approximation of these specifications using a causal dis- 
crete time system 
(3) the realization of the system using finite precision arithme- 
puc 
Note that these three steps are not independent of each other. In 
this thesis we are mainly interested in step 2 and to some extent on 


step 3. 


ιών OF DESIGNING THE DIGITAL FILTER 

An important class of techniques for designing infinite impulse re- 
sponse filters to be realized recursively is based on a transformation 
of a continuous time filter. This class comprises at least three tech- 
niques. 

1. Impulse Invariance Method 

The impulse invariance method, also called standard z transforma- 

tion (or standard z), is a technique in which the impulse response of 
the derived digital filter is identical to the sampled impulse response 
of the continuous time filter. If the impulse response of the filter is 


h(t) then the sampled impulse response will be 
* 
h (t) = h(t)* ó> (t) (20 


where Sr t) is the sampling function and is defined by 


1 t=nT 
Γ n t#nT where n=0,1,2,... 





* 
It can be shown that the Laplace transform of h (t) will be 
E 
& E. 1Κ2π 
H (s) = Д ΗΝ. h(s + 182) (2.2) 
Ka -o 


and the impulse invariance response of the filter will be 


σι ο αι 





7 = е (Zero) 


from the relationship z = cn it 1s seen that the strips of width = in 
и l To map into the entire 2 plane as depicted in Figure 2.1, the 
left half of s plane strip maps into the interior of the unit circle, 
and the imaginary axis of the s plane maps onto the unit circle in such 
a way that each segment of length £T is mapped once around unit circle. 
Thus, the mapping is not a one to one mapping, and hence it can easily 
be shown that the impulse invariance response is only satisfactory if 
H(s) is band limited. And as in most cases if H(s) is not sufficiently 
band limited H(z) is an aliased version of H(s). Therefore the tech- 
nique is used for narrow band filter design or else H(s) is broken into 
cascaded subsections of first order and second order sections with 
esard tılters ın between, which in some cases 1s a tedious job. Also 
it is clear from equation (2.2) that due to the = multiplier, File 
digital filters derived using the impulse invariance method have a gain 
approximately - to that of continuous time filter, which should be taken 
into account in the design. 
2. Matched z Transform 
This procedure is based on mapping the poles and zeros of the 


I 
continuous time filter by the substitution of (5-51) > (l-e 1, 1). 











This means that poles of H(z) will be identical to those obtained by 


impulse invariant method; however the zeros will not correspond. 


3. Designs Based on Numerical Solution or Analogue Differential 


Equations 
a. Conventional Digital Filter Design 
In this method, the differential equation of the analogue 
filter is approximated by a recursive equation, which is the standard 
procedure in Numerical Analysis. There are three basic numerical inte- 
gration techniques, namely, 
e Euler forward integration 
(11) Euler backward integration 
''' linear integration 
All these techniques plus many others are described fully in the litera- 
ture. And there is no need to go into details for all of these tech- 
niques, but because of the importance of bilinear transformation we de- 
scribe it briefly. 


The approach uses the algebraic transformation 





ЕТ 
s = [i2 


2 
з A (2.4) 


to derive the system transfer function of the digital filter such that 


H(z) = H(s) 





2 
= = Ux 5 
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Mis transtormmetion has the effect of mapping the entire left half s 
peto che inside of the unit circle and entire right half of the 

s plane into the outside of the unit circle as shown in Figure 2.2. This 
results in a nonlinear warping of the frequency scale according to the 


relation 


жч 
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DNI. Mapping of f(s) into f(z) as per relationship z = ез1. Note 
that the area between the strips n/T and -m/T map into the 
entire z plane in such a way ne area on the left half s 
plane strip maps into the inside of the unit circle and the 
area on the right half s plane strip maps outside the unit 
circle. The process is repeated infinitely many times thus 
mapping or transformation is not one to one. 

s plane 


1 plare 
imaginary axis 


ES imaginary axis 
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Fig. 2.2. Bilinear mapping of f(s) into f(z) as per relationship s = 
2(z-1)/T(z*1). Note that the transformation is a one to one 
mapping. 
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c NE cd 
m = tan E (2.6) 


where ш. 15 the critical frequency of analogue filter and gts the 
critical frequency of the digital filter. Because of this warping of 
the frequency scale this design technique is most useful in obtaining 
digital design of filters whose frequency response can be divided up 
into a finite number of pass bands and stop bands in which the response 
is essentially constant. Figure 2.3 shows the frequency response of an 
analogue filter and its approximated digital frequency response using 
bilinear transform techniques. 

Brem Prgeure 2.5 it is obvious that if we require the critical fre- 
quency of the digital filter to be say at 0 ¿then we have to frequency 
scale the critical frequency of the analogue filter by a factor 


factor = Bar, 2.7) 
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Typical frequency selective analogue filters are Butterworth, Chebyshev, 
and Elliptic filters. Note that all these filters have closed form 
analogue design formulas and by the use of bilinear transformation we 
can easily get approximate closed form digital filter algorithms. 

A Butterworth filter is a monotonic in the pass band and in the 
stop band. 

A Chebyshev filter has an equiripple characteristic in the pass 
cand and monotonic in the stop band, or vice versa. 

An Elliptic filter is equiripple in both pass band and stop band. 
Clearly these properties will be preserved when the filter 1s mapped to 
a digital filter with the approximated bilinear transformation as shown 


риге 2.5. 
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E frequen 
Analcgue passband Ην 
characteristic 


на) | 
Approximated 
digital pass- 
band charac- 
teristics 


| 


νε 2.5. 


Transformation of analogue filter H(s) into digital filter 
using bilinear transform S - 2(z-1)/T(z«1). 





b. Wave Digital Filter Design 

This technique basically uses the bilinear transformation of the 
analogue to digital design exactly the same as part a iii, but the 
attempt is made only on LC filters having input and output terminating 
resistances of К, and R}. 

In this technique each reactive element in the analogue ladder 
structure is transformed into a two port signal flow structure using the 
bilinear transformation and wave flow techniques, and at the same time 
matching the port one impedance of the derived two port structure to 
the port two impedance of the previously derived element, thus 
eliminating mismatch between the succeeding elements. The details of 
this and its implementation are left for Chapter IV. Note that the 


idea is a very basic one and can be applied on many varieties of common 


filters. 


D. CHAIN OR A B C D MATRIX THEORY 

The analysis of any passive linear network with several inputs and 
outputs can be done in many ways. The most usual ones are signal flow 
graph, system matrix equations, input/output algorithms, transform matrix 
equations, etc. However, for systems of order higher than 2, most of 
the above analysis becomes tedious and prone to mistakes due to system 
complexity. A most useful and convenient method of dealing with a com- 
plex system is, whenever possible, to break the system into several sub- 
systems and interrelate these subsystems by a chain matrix, thus allow- 
ing each subsystem to be analyzed and investigated separately one at a 
time, without even thinking about the rest of the system. To illustrate 
the point consider the network N of Figure 2.4 with wave inputs 8] and 
a_ and wave outputs b. and b.. The relationship between port 2 and port 


2 1 2 
1 for this network can be written as 
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| Network N q 
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System boundaries 


Fig. 2.4. A general two port structure with inputs a, and a, and 
outputs b, and б», into and out of the system boufidaries 





respectively. 
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| | 
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Subnetwork Subnetwork 
N Nə 
| | 
System boundaries 
Fig. 2.5. A general two port structure with inputs a and 


and a 
outputs b, and b, into and out of the system boundaries 
respectively. Nöte that the system of Fig. 2.4 ıs 
equivalent to the system of Fig. 2.5 if the inputs 
(a, and a.) and outputs (b; and 92) are exactly the same. 
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Now if we can break network N into several subnetworks inside the dotted 
line without touching the input and output ports, the resulting network 
will be exactly the same as network N. Note that the networks Ny and N. 
resulting from network N do not necessarily have to have equal subsec- 
tlons — 1n Fig. 2.5. Note also how the consistency in the wave- 


flow direction is maintained. 


iis trom Fig. 2.4 
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Note that equations (2.8), (2.13) are identical. This simple example 
clearly demonstrates how a system when made up of only two simple sub- 
systems N) and N, which are easy to be analyzed each separately, when 
combined into the system N becomes a complex system, and very difficult 
to get analyzed. 

Note that in the first case it is a very easy job to analyze subele- 
ments Ау, BJ; C and D; ОГ system N; or that of the system N}. While 
clearly the analysis of the element A = A r (e cC. Of the system 
N will not be an easy job, and in most cases when the system is made up 
of more than 2 subsystems it is a tedious if not an impossible task. 

Thus equation (2.12) clearly demonstrates the fundamental property of 
the chain matrix analysis. 

Whenever two or more than two pairs are comnected in cascade the chain 
matrix of the composite network will be the product of the individual 
chain matrices. Since in general matrix multiplication is not commutative, 
it is important that the matrices be multiplied in the same order as the 
Circuits are cascaded. It can easily be shown that if all the individual 
Matrices are reciprocal the composite two port will also be reciprocal. 
Б ТОМЕЕРТ OF DELAY FREE FEEDBACK PATH, OR DELAY 

FREE LOOP IN A TWO PORT NETWORK 

As mentioned in section II-D the relationship between port two and 
port one of a subsection of a general causal system such as that of Fig. 


2.5 can best be described by the equation (2.8), i.e, 


Б A By en 

= АЕА) 
a, | а IA δι 
where b), a, are the output and input quantities at port one and b,', 
as! are the output and input quantities at port two. This equation can 


би 





be rewritten in terms of input quantities in the following form 


“wass 
саады a 





bi a] 
= (s) 
D Бо 
or we can write equation (2.15) in the simpler form of 
b; £, (2) f, (z) 8] 
Ξ (2.16) 
be! (2) Аи 
2 | | 4 ες 


Note that for a causal disital system the values E; (z), EXE f, (z) 


and f,(2) are of the fom 


a, * at + --- +а;2 + --- az 
f(z) = - =. (2.17) 
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Note also that when written in terms of positive exponents of z, the 
order of z in the numerator must be equal to or less than the order of the 
denominator for causality. With this in mind the iterative equations 


derived from equation (2.15) are 


b, @) = 0484 (n) + 8,a,' (n) * weighted values of past inputs at 
| port one and port two plus weighted 
values of past outputs at port one 
о) 
апа 


b,'(m)= а.а, (п) + 554, ' (п) + weighted values of past inputs at 


port one and port two plus weighted 
values of past outputs at port two 


(2.19) 
Note that a @, 91) 8, are merely weighting constants, and equations 
(2.18) and (2.19) are both causal and realizable. 
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Now if the input a,' (n) at the port two is a function of b,'(n) of the 


pore two, 1.e€. 


a,' (nj) = Kb,'(n) + weighted values of past input and 
Р ü output values at port two (2.20) 
which is the case for most cascaded two port systems such as the one 
shown in Fig. 2.5 the equations (2.18) and (2.19) further reduce to 
b- (n) - 044 (n) * 8, Kb,' (n) + weighted values of past 


inputs and outputs at 
port one and port two (22 21) 


and 


b,'(n)= a,a, (n) + 3,Kb,'(n) + weighted values of past inputs 
at port one and port two plus the 


weighted values of past out- 
puts at port two (2.22) 


Note that the iterative equation (2.22) is not causal since for the cal- 
Enfation of b,' (n) it requires b,'(n) which is not possible. Thus for 
the causality either K must be equal to zero or 32 mist be equal to zero 
which are two distinct cases. 


Case 1 


By considering Fig. 2.5 it can be seen that a^, Ὀ are merely 
the port one quantities of the second stage. Thus if we are going to 
consider the first stage only we cannot force K to be equal to zero. 
Thus the only variable left is 85 and by making 35 equal to zero we can 
make equation (2.22) a valid equation. This condition is referred to 
the case of no internal delay free path from a,' to b,’ or no delay free 
path in port two. Thus to implement this condition the order of the 


B. 
numerator of the function £ (z) = 12 must be at least one order lower 
D 


1 
than the denominator. 





AS it was noted in the Case 1 since 82) 92) are merely port one 
quantities of the second stage.If we can make 51 ' (n) of the second 
stage to be independent of the a,'(n) of the second stage and only 
dependent on the past inputs and outputs of the port one of the second 
stage, then we have actually managed to make a5' (n) of the first stage 
to be independent of b,' n) of the first stage. This can be done by 
forcing a, = 0 in the second stage. This condition is referred to 

the case of no internal delay free path from a, to b or no delay free 

path in port one. Note that to implement this condition the order of 
Ay (21D, (2) -3, (23€, (2) 

the numerator of the function Е; (2) = = must be at 


least one order lower than that of the denominator. 
F. SUMMARY 


In this chapter we have revlewed briefly the ground work required 
for the matched two port wave digital filter theory and design. 

The contents of this chapter are used throughout this thesis. No 
particular mention of any reference has been made since most of the 
subjects discussed are well established and details can be found in most 
digital filter design handbooks or papers. It is worthwhile to note 
that the concept of delay free feedback path, though important, in most 
papers reviewed were merely stated without any proof. Thus in this 


chapter an attempt was made to prove it in the most general sense. 






JI (Gebers DISCUSSION ON SENSITIVITY THEORY 
ΤΕΕ DIGITAL FILTERS 

A. INTRODUCTION 

The main intent of this chapter is to start with the low sensitivity 
properties of the doubly terminated analogue LC structure and then extend 
this property to the wave digital filter. Later on in the chapter we 
explain briefly the development of wave digital filter theory. It is not 
the intention to go into details, but merely to give an overview of pre- 
vious works for which full development is available in the references. 

mother Objective Of this chapter is to briefly discuss the different 
wave digital filter structures and algorithms which are all called wave 
digital filters, each structure having its own characteristics and limita- 
tions. The newcomer to this field will be astonished by the many differ- 
ent structures and algorithms which are present in the literature; all 
of them are offered under the same name of wave digital filter. 

Later in this chapter a natural development of the wave digital 

filter theory is traced from the initial conjecture of Fettweis up to 


the present day state of the art. 


ОЕ А ТИТОМ OF SENSITIVITY 

The term sensitivity of a certain filter structure has a broad mean- 
ing, and the literature is full of different definitions for sensitivity 
functions meeting a specified requirement. But in general all of the 


definitions end up more or less to the following: 
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As the filter element values or coefficients are varied about their 
true value, the filter pole and zero locations are shifted correspond- 
ingly, causing a change in the magnitude and phase characteristic of the 
filter. A mean shift in the characteristic can be computed on a point by 
point basis and used as a figure of merit, but it may be meaningless 
when determining whether or not the original filter specifications are 
met. Therefore sensitivity must be interpreted in terms of several 
factors to include such parameters as bandwidth, cut-off frequency, 
ripple in the pass and stop band, etc. Several sensitivity functions 
are most commonly used. They are 

(i) Q sensitivity and pole frequency sensitivity 

(ii) Root sensitivity 

(111) Coefficient sensitivity 

(iv) Frequency response sensitivity, etc. 

Note that most of these sensitivity functions are meaningful and for 
most cases there 1s a relationship between most of these sensitivity 
functions with each other [12]. 

Since in this thesis we are mainly interested in the parameter or co- 
efficient sensitivity, thus no discussion of the other sensitivity func- 


tions will be made. 


C. DEFINITION OF PARAMETER SENSITIVITY 


For the function H(C,, Со... [EE ees " where H is a function 


1? 


of multiparameter (Ci, E С.,..., С), the parameter sensitivity 


gru ος 1 n 


H of the system to any change of variables or elements C; is defined in 


Several ways 


C3 


1) S= — 3 m^ logarithmic or (3.1) 
Т factorial sensitivity 
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JR NES = о derivative sensitivity (3.2) 
1 
IDEE ee), semilogarithmic 
Jue H oC - ER EUM 
i 1 sensitivity 


and all of these sensitivity functions are discussed fully in the literature 


021, [15], [14], [15]. 


D. QUANTIZATION ERROR IN DIGITAL FILTERS AND ITS SENSITIVITY EFFECT 

Although tolerance problems which exist in the physical world in the 
case of analogue filters, do not exist as such for digital filters, still 
there 1s an interest in obtaining structures with low sensitivity to para- 
meter variations. There are several reasons for this. The primary reason 
is the fact that the structures with low sensitivity are less affected by 
coefficient truncation. Thus element values or coefficients with 
shorter word length are sufficient for meeting a given specification. 
The second reason stated and proved by Fettweis [10] is that there exists 
a relationship between sensitivity with respect to multiplier variation 
and the round off noise at the output of a digital filter. Thus any 
improvement in the sensitivity would result in a reduction in the corre- 
sponding round-off noise at the output. 

In order to compare the coefficient sensitivity of different digital 

structures Ku and Ng [8] have used a root mean square error criteria in 
the frequency response given by 


N 
ee. = te Y n 2 x m 
es 7 [d 2L, D GO.) [089,17] (3.3) 
1Ξ0 
Note that this criteria is in the frequency domain and is based on the 


deviation of the magnitude of the finite precision output from the ideal 


or infinite precision output. 





Thus ΔΗ(ω:) is defined as 
AH(w,) = [H(o,) - H (w) (3.4) 


where H (w; ) is the magnitude of the ideal output at frequency w> and 
H(w,) is the magnitude of the finite precision output at frequency w;, and 
Ww; J is the weight chosen to reflect the relative importance of error at 


various frequencies,and N is the number of sample points in the frequency 


domain. 

Furthermore in order to equalize the effect of the various coeffi- 
cients, floating point arithmetic is used rather than fixed point 
arithmetic, during the process of rounding off of the coefficient to 


@aemreguired number of the bits. Note that more details and implementa- 
tion of the idea are left for Chapter V where we investigate case studies 


based on bit quantization error. 


E. SENSITIVITY OF LC LADDER STRUCTURES A 

It is a well known fact that the doubly terminated analogue LC 
ladder structures are relatively insensitive to element value changes. 
To prove this fact most researchers have used the theory of the total 
conservation of input and output power quantities in a reactive lossless 
network, the discussion of which is interesting but it not τη τα (Cope 
of this thesis. A very good review is made in this respect by Renner 
and Gupta [7]. 

Thus, it is a reasonable assumption that the digital filter derived 
directly from the topology of a resistively terminated analogue lossless 
structure would have the same favorable sensitivity properties as its 
analogue counterpart. 

To show the low sensitivity of wave digital filters by the conven- 


tional n port scattering matrix network theory, Fettweis [10] defines the 
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instantaneous "pseudopower"' transmitted through a wave port k with port 


input resistance R, by the means of the equation 
(3-5) 


where p, is the total instantaneous pseudopower input through port k and 


a, is the total input wave at port k and b, is the total reflected output 


wave at port k. Similarly he defines steady state pseudopower by the 


relationship 


) (3-6) 


where P. and Ay and B. are the steady state values of Рр аг, b, ВЕРЕ 
tively at an arbitrary frequency "f". In this way he then shows that for 
all wave digital building blocks which are derived from the LC doubly 

terminated ladder structures, the total sum of instantaneous pseudopower 


Mesorbed bv all ports is equal to zero, i.e. 


Yo 


k=1 
where n is the total number of the ports of the wave digital subelement 
or building block and P. is the instantaneous pseudopower at port k. 
Similarly he also proves that the total steady state pseudopower 
absorbed by all building blocks of a terminated wave digital filter at an 


arbitrary complex frequency "f'" is equal to zero, i.e. 
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where Ρις 15 the steady state pseudopower at an arbitrary complex frequency 
f, and n is the total number of the ports of the wave digital building 
block. The low sensitivity property then follows from the zero pseudopower 
relationships. Although power has physical meaning in analogue circuits its 


meaning in digital filter algorithms is nebulous. 


F. PREVIOUS WORKS ON THE DEVELOPMENT OF THE WAVE DIGITAL FILTERS 

In the previous section we discussed the favorable low sensitivity to 
Внтамесег variation properties of LC ladder structures. It was this fact 
that led Fettweis to propose the idea of wave digital filters in 1970 [5], 
derived from the doubly resistive terminated LC ladder structures. Since 
then numerous papers relating to the wave digital filters have been presented. 

To understand the wave digital filters and the state of art, it is 
necessary to summarize the state of evolution of the wave digital filters. 

Partial credit of the development of the wave digital filters 
can also be given to Richards [1], Kuroda, Ozaki and Ishii[2] and numerous 
other researchers who contributed to the development and synthesis of the 
Strip line filters and also resistor-transmission line circuits. 

The basic idea behind wave digital filter development from the beginning 
was to keep the desired low sensitivity of lumped LC resistively terminated 
filters in the already developed analogue domain and transfer it to digital 
domain with very minor modifications. 

Ideal strip line filter and transmission line circuits having inherent 
LC structure is an ideal starting point. To avoid the loading effect of 
@me clement of the LC strip line upon the succeeding and preceding elements, 
Ozaki and Ishii [2] introduced the idea of inserting unity element, or the 


lossless transmission line acting as an ideal transformer (which was already 


developed by Kuroda and known as Kuroda's identity), between the stages as 
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a lossless matching element. 

In a more detailed paper Fettweis [4], in order to arrive at realiz- 
able signal flow diagram, used wave quantities instead of the usual 
current and voltage input quantities to the conventional filter. Actually 
in using the wave quantities as input and output signal, Fettweis made 
use of the already developed scattering parameter matrix theory of the n 
port networks. In order to avoid the loading effect of one section upon 
the other, and thus to eliminate the unwanted mismatch reflections between 
sections, he employed unit element or impedance transformer, known as 
Kuroda's identity [2]. 

Although with the introduction of the unity element a desirable re- 
fete Was achieved in realizing the LC filter structures, this introduc- 
tion further increased the number of required operations (mainly multipli- 
G OS) over the conventional digital filter. 

In order to avoid the unity element several attempts were made by 
different researchers, Sedlmeyer and Fettweis [5], Van Haften and Chirlain 
[6], and others, to eliminate the need for unity elements or impedance 
transformers between the cascading sections. One of the attempts was made 


by Fettweis [5] and later on the method was renamed by Ku and Ng [8] as 


‘the "new Fettweis method." The new Fettweis method forced one of 
the scattering multiplier coefficients in the three port network to be 
equal to 1, thus a two port element was made from the three port element. 
By starting at one port and progressing forward towards the other port, 
and by suitable choice of input impedance of the next section (1.6. 
matching the output impedance of cne section to the succeeding section), 
the need for unity element or impedance transformer matching was 


eliminated. 
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Van Haften and Chirlain [6] next came up with the idea of cascading the 
wave digital three port section directly without the unity element but 
taking into account the approximate attenuation of one section upon the 
next section, thus eliminating the need for unity element and extra 
multiplications associated with it. At the same time this introduced 


errors associated with the calculation of the approximate attenuation of 


one section upon the other. 

S. Erfani and B. Peikari [9] and N. S. Swamy and K. S. Thyagarajan 
[11] at the same time came up with a new approach to eliminate the unity 
element. This technique is the basis of the further investigation in this 


thesis, and details are given in Chapter IV and other chapters. 





ου 


ED. 
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IV. GENERAL THEORY OF THE WAVE DIGITAL FILTER 


A. INTRODUCTION 

The intent of this chapter is to derive the generalized algorithm 
for a two port wave digital filter on a step by step basis using the 
principles of the circuit theory, matrix algebra, and scattering matrix 
wave theory. 

The resulting algorithms are summarized in two tables in section E. 
Also, two 11lustrative examples are given, example one being a simple 
element, suitable for matched source; and example two being a complex 
element, suitable for a matched load. 

It must be emphasized that the techniques used to derive the wave 
digital algorithms in this chapter are more general than the previous 
works, in the sense that only one algorithm is derived for both series 
and shunt element. Also the effect of sampling time is introduced into 
the algorithms. 

B. DERIVATION OF THE GENERALIZED TRANSFER RATIO 

FOR A TWO PORT WAVE DIGITAL FILTER ALGORITHM 

Consider the two port network N, of Figure 4.1b whose inputs are a, 
and a, and outputs are b; and D>. In order to use the chain matrix theory 
developed earlier (Chapter II, Section B-D) we need to find port two 
input and output waves 82292 in terms of port one input and output waves 


a} >b} or vice versa. 
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Eo. 4.1. 


L 


Two generalized representations of a two port network. 

Fig. 4.la representation in terms of voltage and current 

suu. ασ, 4 1ο representation in terms of scatter- 

ing matrix wave quantities. Note that R, is the input 

impedance of the network and R, is the ouitput impedance 

of the network. R, is the source impedance and R, is the 
S L 

load impedance. 


T 





Thus considering the scattering matrix model of Figure 4.1b, we have 


1 T I 
= ΠΠ; 
b 1 -R || 
ш ei – E 
and 
а, 1 К, V. 
= (А6. 
[Usi ΕΗ 


where R) and R, are port impedances of port one and port two respectively. 


Note that a, ,D and a,,b, are voltage waves. Now in order to find 


| 
82202 in terms of 81,51 we have the chain matrix of 


Va A B V. 
= 5) 

I E D I 

lo 2 | 

Beuations (4.1) and (4.3) lead to 
} ү 

а, E К, А B 2 

y! |" Ἢ ἐξα 
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OF 


























a, А R C Iur 
b; A zi RC B = 
But from equation (4.2) 
-1 
ne — 
V 1 R 
2 2 
1 “Ry | 
or => ЕЕ 
V. 5 > 
I, 3282 E К. 
Thus equations (4.4) and (4.5) lead to 
EU Ae A 
A R.C a +R C 
2 2R 2 
A- ΞΡ 
λ RC 1 R C 
D 2р 2 2 
Or 
ај u À 
b; V < 
where 
A+R,C 
| | 
HM + 5 (B+R,D) 
2 2R2 l 
A+R С 
cil il 
À = 5 - 2R> (B+R,D) 
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R D 


x 
(АА) 
RD ДЕ 
2 
Do 
= 
2 
(4.5) 
52 
1 
(4. 6) 
1 | 
үф 
b2 
(4. 7) 
Ho 
(2:735 
(4. 7b) 





У = — ik (B-R,D) (4.7с) 
A-R4C 1 
к= л = R, (B-R.D) (4.7d) 
Thus from equation (4.7) we have 
s (ius =“ El (4.8) 
1 | u 2 
o À 


and the appropriate wave tlow diagram is drawn in Figure 4.2. 


Note that 
V =V -IR (4. 10) 


(4.11) 


Thus equations (4.2) and (4.11) lead to 


22 = Y, + El, 


ο. V. š R51, 


2 
V =- LR 
ES that a, = bod (4.12) 
2ү, 
05 = 155 (4.15) 
where 
NU > 
Ec (4.153) 
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; Ry 
2 T 
Ro >51 Ц u Re ->R 


Network № 


Fig. 4.2. Wave flow diagram representing network N, of Fig. 4.1. 


E 





and also equations (4.1) and (4.10) lead to 


τ τ; 
T 
ПА ° LRS 
so that a, + ӨЫ; = НУ. (4.14) 
where 
5 = Е - 3 (4.15) 
1 S 


Where V. is the input voltage. Thus from equations (4.13) and (4.14) we have 


V po íz) 
2 (1+5) (+ 2 
H(z) - peu nme (4.16) 


Note that for 9 = 0, from equation (4.15), Ry = Ro: 





ME > ox 
and (z) 
DAZ 
OA, 2 


For 4 » 0, from equation (4.13a), К, = Rr and from equation (4.12) a, = 0. 


Thus from equation (4.16) 





be) 
Λη 2 
H(z) = = == (4.18) 
2 a. (z)*9b4 (z) 
pr 
be) 
H(z) = 1*9 2 (4.19) 


2(1 + 9 2 а (2) 


Сл 
ΓΩ 





The general realization with the aid of equations (4.11), (4.14) and 


(4.17) is shown in Figure 4.3. 


C. REALIZATION OF TWO PORT WAVE DIGITAL STRUCTURE 

Theoretically LC passive structures such as ladder, lattice and most 
of the other symmetrical structures can be reduced to a form of Figure 4.4 
which is the "T" form with positive or negative element values. Note that 
NA BC D*matrix,the direction if 1} 1S in the direction as shown in 
Figure 4.4 and 4, £, or 4, can be positive or negative capacitors, inductors 
and/or tank circuits with positive or negative element values. 


The A BC D matrix of Figure 4.4 as per equation (4.3) is found from 


the equations 








E. 2,2, 
V 1 + = Z + 2, + 5| [V 
l 2 1 3 72 2 
= (4. 20) 
z 
© 1 D Ee ЈЕ 
7 52 52 
Thus 
А = 1 + ED (4. 20a) 
22 
Iu a 5 
B an eu + = + 2 (1 20b) 
ee (4. 20c) 
=) 
E 
Ши = = (4. 20d) 
22 





Network Ny 





δι ἃς К-К) 


$ = —— 
КК. Ry +R 





Fig. 4.3. General realization of a two port network N, with input 
and output waves. Note that a. Ten ee input voltage wave. 





E. A general two port "T" network with camplex elements 
21» 255 Z-.. Note that R, is the input impedance of Port l, 
and К Ts е input impedance Or Portes 





Note that for a single series element Z, as per Figure 4.5 and equation (4.20) 


2.20 
ә 
7, = = 
1 2| | V; 
0 i MAT 





and for a single shunt element Z 


Bp 
~ 0 
MT 1 0 у, 
1 
I — 1 í 
T 72 y 
b. m аы 4 


= 
Now we can easily find the values of v, à, u, K. of equation (4.7). 


Thus equations (4.7a) and (4.20) lead to 


A 72) 
+ R. + Z c qe ee AC Io 1 2 3 

1 2 1 5 22 

и = о C 00 (4.21) 
ZR, 


and equations (4.7b) and (4.20) lead to 


ο ο ο ο ο σα 


(4. 22) 





| | 
V —- Ry | |R <— 2 
| ро | пе 
ο occ p e 
ë TEMA ~F P 


Fig. 4.5. ^A general two port with series element only. 
Note that К, апа z are input impedances of 
Port 1 and lPort 2 respectively. 





Fig. 4.6. A general two port with shunt element only. 
Note that R, and R, are input impedances of 
Port l and port 2 respectively. 





equations (4.7c) and (4.20) lead to 


=p. 


τ᾽ ᾱ7 


Ry i 24) (R^ B Zz) 


5 2 
у = 282 (aes) 
equations (4.7d) and (4.20) lead to 
(Z, - Ку) (К - 25) 
Ву + В, - 7 MNA Dx 
Z 1 5 22 
к = ЗК, (4.24) 
Thus the elements of Figure 4.2 will be 
E ^ (75) 
u ENDE I) 
R. M EST a Ry 1 ( Z 5 
2 1 5 2, 
(- * Z4)(R, + 75) 
- +В, + 1 t s *+ — 
~ = (4.26) 
H CR, + 2)) (R, + 25) 
Ir + 7_ + 
1 Жо ES 22 
AV 267 
EN C e d 
Е R Es 
ee, e > 
1 2 j| 3 2, 
О) 
MSN m Toc > 
τ us (4.28) 
+ 7 E 
μιν... е. 
ES — 


Note that for the single series element of Figure 


Ш), (4.27), (4.28) reduce to 
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4 


.5 the equations (4-25), 





TES 


c 


and for a single shunt element of Figure 4.6 the aforesaid equations will be 


== 


L 


ς - № 
u 
_A 
u 
where 
σι 
and 
е7 
i 


nu y 


ο aS 


aj 
e | 


U 1 
e 





τ 7 
ке ai 
ES + R^ ο 
Ry = К, + 7 
^К, 
R, +R, + 4) 
К - К +21 
+ R, t ¿J 


2G, 


> 67 +6, 


26, 


о 


1 


+ G) - G2 


S 


(4. 29) 


(4 .30) 


(4.31) 


(4.32) 


5, 


(4-34) 


(4.35) 


(4. 36) 








D. DERIVATION OF GENERALIZED ALGORITHM WITH NO DELAY FREE PATH IN 
PeRteONG OR SIMILARLY WITH NO DELAY FREE PATH IN PORT TWO 


The wave flow diagram of Figure 4.2 is not suitable to be used in 
the chain matrix equation,since most probably it has delay free paths 
between a, to b; and a, to b, which we have not analyzed as yet,and if so 
enot desirable. Thus by a careful choice of either Ry or 
R, we can force the wave flow from either а, to b, to have no delay free 
Bath, or from a, to b. to have no delay free path; but not both of them 
since we have only one variable to adjust, and that is either Ry Or 
R5. 

The delay free wave flow for the case of no delay free path from 81 
το δι for a three section filter is shown in Figure 4.7 and the delay 
free wave flow for the case of no delay free path from a, to b. for the 
Same section filter is shown in Figure 4.8. 

Note that a close examination of Figure 4.7 and Figure 4.8 still 
reveals a delay free feedback loop at the terminating points. Thus to 
make the sections of practical value, $ must be forced equal to zero for 
Figure 4.7 implying that 

A 


o = Rp E = 0 or Кү = К, (4, 57) 


and for the Figure 4.8, 6 must be forced equal to zero implying that 


А 
1 S In x 
mm") or К, = (4.58) 
S 


Thus the updated and useful delay free wave flow diagrams would be as per 
Euupres 4-9 and 4.10. Note that in the first case we have to start at 


load end and work our way through to the source end, and in the second 
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case, l.e. Fig. 4.10, we have to start at the source and work through to 
the load end. Alternately we can start from both ends and progress 
towards each other. Where they meet we have to introduce a transformer 
or unit element such as Kuroda's identity, but this is not really 


recommended since it introduces further complications. This approach 


is not desirable and is to be avoided. 

КО TI DC aid of formulas (4.25), (4.26), (4.27), (4.28), 
we can cascade as many elements in a section as complex as Figure 4.4 
or as simple as Figure 4.5 or Figure 4.6. 


Note that the element values 21» 22 and 3 mist be in s domain, 


and tor digitization we use the bilinear transformation 





(л 
II 
-ἠ[το 


Е (4.39) 


where T is the sampling period of the digital filter. 


Example 4-1 
Series L, no delay free path from a, to By еле ок азирге 
element : 


With reference to Figure 4.2, v/u must have no delay free path, and 


fae aim is to find RI given R,, the load end terminating resistance. 


DES 
From equation (4.50) 





v ĉl Sige? 

р 2 + RL +R, 
and 

п-т | 
Tnus 


2b a _ . 2L. 
E te TV 
DENM ο ο. (4.40) 
2L 2L 
от 
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m to be with no delay free path from equation (4.40) 


2L 


pes 9 
hich leads t R, = R nc 
which leads to 1 = № + 7 (4.41) 


Brom equations (4.29), (4.30), (4.31), (4.32), (4.38), (4.39),(4.41) 


we have 


(4.42) 


= | 
| 


- σ(ζει) 
(2+0) 





(4.43) 


ο а a (4. 44) 


λ ] 
a r Со) 


IN] 


Z(1-c) 
Ve 


Thus the output equations for series L with no delay free path from a, to b 








1 
аге 
_ o-l zei 
DE 5.7 
b = geb) a E 2 (1-9) а 
2 z*g 1 240 2 (4.46) 
a -1 1 
1 + 02-1 1+с2-1 
-1 

b. - σ(1:-7 a, р л. a, (4.47) 

l + oz l*cz 





with G = 


prp 


and the iterative equations, corresponding to (4.46), (4.47) are 


b, (n) (c-1)a, (n-1) * а, (п) + a, (n-1) = ob, (n-1) (4.48) 
b7(n) = ca} (n) + ca, (n-1) + п. - ob, (n-1) (4.49) 


With zero initial conditions, i.e. 


X) = 0 
Х = 0 
о 
X, 7 0 


The equations (4.48), (4.49) can be rewritten as 


b; = σ(χη - xs) a х, + a, (1750) 


b2 = o(a, - à, + Xx) - x4)+ a, (4.51) 


and updated equations 


EE ST (4.52) 
хә = by (4.55) 
Xz = a, (4.54) 
x4 = b, (4.55) 


Note that there is only two multiplications for each iteration. 
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Example 4 > E 


ing the composite series and shunt element with 


An example of design g mp 


no delay free path on port two. 


As an example of a composite series and shunt algorithm, let's con- 


Sider the two element section of Fig. 4.11. The aim is to derive the 


wave flow algorithm with no delay free path from a, to 55. From equation 


(4.28) we have 


і 
т 


where Ay 


Cy 


With reference to Fig. 4. 


have no delay free path, 


Ry - К, + SL - (R; ου R, sC 























Кр + К, + 51 + (Ку + 51) К 5С 
2 
2 
КЕ s(L + R, RC) ES (R4LC) TE 
5" T GT 
7 
uU ym 
FERE 2 
= EE E 
τ: 
_ | E 2R R.C | 4R,LC 
PES ee т z 
AR LC 
1 
2L 281826 4R,LC 
=R -R - — + - = 
р T т? 
A Bi 2R R.C | 4R4LC 
| 2 T T T? 
τι 
πι... eo 
y j 2 т? 
_ m. Au ZR Е.С | 4R4LC 
т 72 
2 for no delay free path on port two, - must 


thus 


run 


— 
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| 
port 1 ---------- --------- porto 


Figure 4.11. Composite series and shunt element of example 4.2. 
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This leads to 





Thus from equations (4.57) and (4.28) we have 


























у 8,27. + pz 
s E =) 
er 817 + yı2 
where 
RC RE 81,^С 
ο а, 
Т 72 Т τ 
52 < NOE on (5 
Ἢ ΓΤ n 2 
2 
2R C 1. qs 
" T Т 13 
2 = RC 
< 21, П τς 
Ry + PIC + a E 
R ^C 2 
T SL^C 
2R + == * : 
- T | το 
i Е REG 
2L I“ AC 
(R, + >=) T m 
Ry 
ΠΤ 2R.C 
2L J ni 
E 7 


and from equations (4.57) and (4.25) we have 








-1.2 

= 

u 
J + 812 1. Y? 4 

where 

д | 

: RC 41 
1 + T Mr 





(4.57) 


(4.58) 


(4. 59) 





and from equations (4.57) and (4.27) we have 


AV 84 (129 
rac Кс озен ο (4.60) 
T 812 + у12 


where 
Ri 
G R 
p= 
ΗΠ 
and from equations (4.57) and (4.26) we have 


А ао. (4. 61) 


Thus the wave flow algorithms for the composite element with no delay 


free path from a, to b. Will be 


ey? о + acne 
65 = да. + oo a (4.62) 
: -1 -2 -1 - 2 2 
px 812 A Tor 8, 2 + ¥42 
-B52 - Y> a, Q«z 7 D* 
И I, d (4.63) 
1 E -2 1 E E 
D 592 + yı2 1 8,2 + Y 2 
or 
= 7 - E — - ‚64 
b, (n) 8а. (п) + 28-а (п-1) в-а (д ) + За. (п 1) (4.64) 


* Y5985(n-2) -B4b5(n-1) - y4b5(n-2) 


val 





b, (n) m -8584 (ПЕ) -Y584 (n) t 8 ,a, n) + 28 a, (n-1) + 8 ,a,(n-2) 
-В161 (n-1) -y, b, (n-2) (4.65) 
and the iterative equation corresponding to equations (4.62), (4.63) with 
113] conditions set to zero, i.e. 
х. = 0 1=1,8 


will lead to the iterative equations of 


= 22 > Е 
π.Χ. - ΒΚ - κ eo 
by = -Y28 - 82%] * 84(a5 + 2x3 + ху) - BiXgs - үүх (4. 67) 

with updated equations of 
үа 
οσα 
эте; 
Xy = x; (4. 68) 
Хе т b; 
us 
Xz = bo 
2398 727 
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Note that there are 12 multiplications with six coefficients, i.e. a's, 
B's, y's. It is interesting to note that since in this example we 
originally had two elements, 1.e. Series L, and shunt C with input in- 
pedance Ву. Thus we expect the a's, B's, and y's to be functions of Ry» 
L, and C only. This can be easily shown by identifying two new variables 


a and 8 such that 





2L 
ir T 
= ^L 
NT 
1 
8 = 
йаш ше 
jt T4 
Then 
=, 2 (8) 
Lucem c Кт 


51 ιν ζα + Bla+l) 
Ү1 = gd -a) 
32 = ] - Salzen 


Y2 о 8(a-1) 


B4 = l-a 


5 i 
R> = = is more general even for the сазе Eg zi =( which makes a=1 and 


l-a=0 and if we use R, = Eu. К, becomes which is indefinite, while 
К, = a 1S not indefinite. 
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By inserting these new values into the iteration equations (45.66) and (4.67) 
the number of multiplications reduces from 12 to 9 which is significant. 


Thus the new iterative equations from (4.66), (4.67) will be 


b; 7 85-84 XQ *2x- ЈА x à (-a5-x4*2(84-X.*Xc)) (4.69) 


+ 


8 (a4 *X1-Xc-Xg) + aB(-a1*X4-Xg*Xe) 


о" 
N 
| 


= x3*Xy -X4*2a(X--X,) т 8 (a4 *2X4 *X-X5-X4-X- -X4) ΠΠ 


Xx 


αβ (x: -Х7+Х9) 


Note that the multiplication a8 and 2a does not enter into the iteration 
and they can be premultiplied before the starting of the iteration process. 
Notation 
For latter analysis in Chapter Y we call these latter set of equations 
I.e. equations (4.69) and (4.70) as "reduced parameter complex wave digital 
algorithms" and the original set, i.e. equations (4.66) and (4.67) as 
"complex wave digital algorithms". 
Although theoretically it is always possible to find new variables 
Samal to the mumber of original L's and C's, in practice it becomes a 


tedious job for more than two variables, and it is not recommended. 


E. TABULATION OF SIMPLE L AND C ELEMENTS IN SERIES AND SHUNT 

In a similar way as examples 1 and 2, the wave flow equations for 
the single L and C elements for both cases of series and shunt are 
tabulated in Tables 1.1a,b and +4+.2a,b. 

Note that Tables 4.la, 4.1b are for the case of no delay free path 
from a, to δη» and Tables 4.2a, 4.2b are for the case of no delay free 


it 
path from a, to b.: 





Table 4.la. Algorithms for simple L and C elements with no delay free 
path from a, to b, with sampling time T. 
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Table 4. 1b. Iterative algorithms for simple L and C elements with no 
delay free path from a, to b, with sampling time T and 


with initial conditions i.e. ХХХ. > set to zero 


Element ЕСТ а С Прота Ст Remarks Updated Values 
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Table 4.2a. Algorithms for simple L and C elements with no delay free 


path from a, to b, with sampling time T 





Element Algorithm | Remarks 


ο] 


Series C 
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Table 4.2b. Iterative algorithms for simple L and C elements with no 
delay free path from a, to > with sampling time T and 


with the initial conditions I.e. X;,X5,X.,X, Set to zero 


Element Iterative Algorithm | Remarks Updated Values | 


| 
| 
| b 
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F. FURTHER INVESTIGATION OF THE ALGORITHMS WITH NO DELAY FREE PATH IN 
PORT ONE, WITH THE ALGORITHMS WITH NO DELAY FREE PATH IN PORT TWO 


A closer look to the two kinds of algorithms derived, i.e. algorithms 
with no delay free path from a) to b, and the algorithms with no delay 
free path from a, to b,, will reveal that both algorithms would behave 
very much like their equivalent LC structure matched to load or matched 
to the source as per the reciprocity theorem. Note that this can also be 


proven by writing the chain matrix equations 








E E a 
= For no delay free 
b, |" κ 7 path from a, to b; 
Or 1 Ky νι 
ENS A al 
or ja oe 
2 т 1 1 
E With no delay free from a, to b; 
E A у 
So if we interchange ay and 82) and also by and b, we would get 
a oli 
1 1 
E With no delay free from a, to b, 
bj a 
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Note that if we calculated the chain matrix for no delay free path from 


a, to b, i.e. the matrix 


a. HJ λα 5, 
b; У, к) а, 
Then the ti -1 
en vO matrices K1 νι н, A, 
and 
Az H] У, А 
would be identical, i.e. 
el an Ἡ2 Az 
Em ΠῚ 
E € 
— Y — .P Q ν 
ΠΠ πι ΠΠ 2 2 
81 
Thus Mose Е 
ЕУ 
À, = 2 
E Ч 
2 A191 x: Satan 
K `1 
2 EI ~ ληνη 


G. DESIGN EXAMPLE 
As an application of the illustration of the techniques developed, 


a Seventh order low pass Chebechev filter with the given specification in 
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Appendix lÀ was designed and implemented in both time domain and fre- 
quency domain. As a comparison the same filter was designed using the 
conventional digital techniques. The computer program and output results 
for both cases are also given in the same appendix. As can be seen, they 
are completely identical in both time domain and frequency domain, thus 
assuring the validity of the theory developed in this chapter. 

Note that in this example the double precision arithmetic was used, 
thus ensuring infinite precision for multiplier coefficients, to come up 


with exactly the same answer. 


H. SUMMARY 

In this chapter we developed the generalized algorithm for a two port 
wave digital filter step by step from the first principles of circuit 
theory and background given in Chapters II and III. The generalized LC 
system considered was an LC "T" Section. For a m section the procedure 
would have followed in the same way with the exactly Similar end results. 

The delay free path does play an essential role in the algorithms. 
Thus emphasis was given on this matter and whenever possible explanatory 
delay free paths were illustrated in the diagrams. To make things as 
clear as possible two worked examples are given, one with a simple 
section with no delay free path at the terminating source port and 
second example with a simple complex section with no delay free path 
at the terminating load port. Then we investigated the relationship 
between matched source algorithm with matched load algorithm of a given 
system. Finally a practical design of a wave digital filter using the 
already derived algorithms were implemented and the complete results 
and computer program are given in the Appendix 1В,-Е. 

Note that no reference is given at the end of this chapter. The 


only relevant references would be references [9], [11] of Chapter III. 
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V. SENSITIVITY ANALYSIS OF THE WAVE DIGITAL FILTERS 
DUE TO TRUNCATION IN THE NUMBER OF BITS 
A. INTRODUCTION 

The intent of this chapter is to investigate the behavior of wave 
digital filters to quantization in the number of bits in the multiplier 
coefficients. The wave digital filter has been credited in most of 
the published literature as being minimum sensitive to multiplier co- 
efficient variations, and thus conjectured to be optimal in requiring 
the fewest number of significant places to achieve a given specifica- 
Шоп. 

In this chapter, for a given analogue resistively terminated LC 
filter, three different wave digital filter algorithms are computed, 
and their behavior with respect to truncation in the number of bits 
of the multiplier coefficients are investigated. The results are com- 
pared with conventional digital filter design. It is important to note 
that for all the algorithms under investigation the bilinear transfor- 
mation is used to make the comparison meaningful. 

In order to avoid any confusion in this chapter we use the term 
"multiplier coefficient" for the coefficients of wave digital filter 
algorithms and the tem "polynomial coefficient" to denote the coeffi- 
clients à5» 85,777, and бе» бе 84 of H(s) or H(z), of a given 


analogue or digital filter, respectively; where H(s) or H(z) are of 


the form 
1 
H(s) = > 
s +a s +a ss +a s tas +a,s”+a, sta, 
H(z) = —— T 21) 
Пор ρου 
2 +562 +657 +512 +052 +552 +61 2+0] 
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It is also important at this stage to note that since wave digital 
Piecer multipliers are derived directly from the original L's and C's, 
they can be thought of as a modified or matched L or C elements. Since 
wave digital multipliers are not derived from the polynomial coeffi- 
cients ofthe transfer function of filter but from L's and C's of the 
analogue circuit,it is meaningful to compare the effect of bit quanti- 
zation error on wave digital filter multipliers with the effect of bit 
quantization error on L's or C's in a conventional digital filter 
algorithm, and not to compare with the effect of bit quantization error 
on the filter polynomial coefficients. Based upon this argument, the 
original filter L's and C's are quantized in the conventional digital 
filter algorithms and used for comparison purposes. 

B. COMPARATIVE STUDY OF THE ERROR DUE TO QUANTIZATION IN 

THE NUMBER OF BITS IN VARIOUS DIGITAL FILTER ALGORITHMS 

Coefficient errors introduce perturbations in the zeros and poles 
of the transfer functions, which in turn manifest themselves as errors 
in the frequency response. In Chapter III Section B we defined general 
sensitivity, and in Section C we discussed parameter sensitivity, and 
in Section D sensitivity due to quantization in the number of bits. 

The comparison criteria we have used in this chapter is the root mean 
Square error criteria given in equation (5.5) and repeated here in 


greater detail, i.e. 


N 


E 24% 
E y A SCA 1 | 
1=0 


where W(o;) is the frequency weighting function. 


ΔΗ(ω:) is the difference in magnitude between infinite 
precision algorithm and the algorithm obtained 


in using a finite number of floating bits in the 
filter multiplier coefficients or component values. 


N is the number of sampling points in the frequency domain. 
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ως 15 the frequency associated with the sampling point i 


It should be pointed out that in floating point calculations, the man- 


tissa of the floating number in the binary system is truncated to the 


specified length using sign magnitude representation. 


However the word 


length requirements of sign and exponent are not counted, since they do 


not enter into the calculations. 


in total nine different normalized low pass seventh order filter 


algorithms were chosen for investigation from the handbook of filter 


synthesis [1]. 


1) 
ii) 
iii) 
iv) 
v) 
vi) 
vii) 
viii) 


Ex) 


These filters are 


Seventh order 


Seventh 
seventh 
seventh 
seventh 
seventh 
Seventh 
seventh 


seventh 


mach filter was 


a) 


b) 


c) 


d) 


.9 db ripple Chebyshev low pass filter with ' υπ 


order .1 db ripple Chebyshev low pass filter with RL, κ 


order Butterworth low pass filter with R.-1, К =1 


order 
order 
order 
order 
order 


order 


.9 db ripple Chebyshev low pass filter with R.-0, 
.l db ripple Chebyshev low pass filter with R.-0, 
Butterworth low pass filter with R_=0, R, =1 

.5 db ripple Chebyshev low pass filter with R -10, 
.l db ripple Chebyshev low pass filter with R.-10, 


Butterworth low pass filter with R.=10, R71 


designed using five different algorithms. 
> 


wave digital filter with simple sections, matched to the 


load (i.e. with no delay free path in port two) 


wave digital filter with complex sections, matched to the 


load (i.e. with no delay free path in port two) 


wave digital filter with complex sections but with reduced 
parameters and matched to the load (i.e. with no delay free 
path in port two 


conventional direct digital filter of the form 


84 


Қ =1 


R, 51 


R =l 


Rr =1 





HZ) = 
; | l+a 271 + а ^2 т - 5 -4 -5 -6 
1 2 4-2 -7 


(e) conventional cascaded digital filter of the form 


mer 
Pez) = ra) 


Е. 
5» = 0 = 
(1+a,z ) (1+0.,2 + 8,7 ) (1*azz I + 8-2 2) (1+a4z 1 + 8,2 7) 


Note that all the filters are low pass and have zeros at z = 1. Thus 
there is no necessity for pole/zero pairing in order to optimize the 
response of the cascaded conventional digital filter. Thus all these 
examples are in a naturally optimum arrangement. 

Double precision arithmetic used throughout in order to get 
practically the same result for all five algorithms with infinite pre- 
cision arithmetic for the coefficient or component values. 

The root mean square error in the frequency response for various 
fllters under examination were found, and tabulated in Tables 5.1-9. 
The corresponding graphs are drawn in Figs. 5.1-9. Note that in order 
to avoid congestion the error graphs of conventional cascaded digital 
filter are not drawn. 

A sample computer program for each of the five algorithms used is 


given in Appendix 2-A,F. 


C. STUDY OF THE RESULTS OBTAINED 
The study of the graphs for all different algorithms clearly indi- 


ве the following points: 
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(1) The logarithm of rms error versus the number of bits £or the 
direct conventional digital algorithms follows a straight line for 
all seventh order filters and they all pass close to the point of 
error - l at bits - 0 with a slope of approximately 3 db per bit. 

Also the rms error of the conventional direct digital filter is never 
worse than the other filter algorithms. 

(2) For filters with source resistance К, = 1 and R. = 0, the rms 
error graphs for the simple wave digital filters are exactly the same 
as that of direct conventional digital filter, with the same slope. 
Ber filters with Ro = 10, the rms error graph for the wave digital 
filter is worse than the conventional direct digital filter. This error 
behavior is very interesting and shows that wave digital filters do 
have a higher sensitivity to filter terminating resistances. 

(3) The rms error graph for the complex wave digital filter for 
ἃς = 1 and ἃς = 0 1s worse than that of simple wave digital filters, 
With approximately the same slope. 

(4) The rms error graph for reduced parameter complex wave digital 
filters with п апа ЕЕ 0 in some cases is slightly worse than 
that of complex wave digital filter, but with approximately the same 
slope. 

(5) The rms error performance for all wave digital filters with 
R. = 10 is approximately the same, after first four or five bits, and 
Ege rms error graph tends to follow a straight line. Departure from 
the straight line is more noticeable for the simple wave digital 
filter than the complex wave digital filter or the reduced parameter 


complex wave digital filter for the first six or seven bits. 
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(6) Although, ideally we expect the two conventional digital 
filters (i.e. the direct and cascaded digital filters) to have the 
same rms error, the factorization process used for the denominator 
introduces some slight error in the pole locations, which in some 
cases changes the rms error value. 

D. COMPARISON OF RESULTS WITH EACH OTHER AND WITH 

THOSE PUBLISHED IN THE LITERATURE 

The results obtained in this chapter and in Section C indicate 
the following points when compared with each other and with results 
published in literature: 

(1) For all purposes the digital filter derived directly from 
analogue doubly terminated LC structure is the least sensitive struc- 
ture with respect to filter component values. This fact should not 
be confused with the sensitivity of the digital filter to polynomial 
coefficients. The foregoing also apply to band pass or band stop 


filters derived from analogue low pass LC realization. 


(2) It is worthwhile to consider the structure of wave digital 
filters compared with the structure of analogue LC doubly terminated 
fllters for cases under study (i.e. seventh order low-pass filters). 
Ihe total number resonances in the analogue LC filter shown in Fig. 
3.10,is ten. These are evenly distributed resonances and can be 
compared with the total number of delay free feedback paths in the 
wave digital filter, which for the case of simple wave digital filter 
is eight as shown in Fig. 5.11. Note that the delay free feedbacks 
are localized towards the source for wave digital filters with no delay 
free path in port two. Also, as shown in Fig. 5.12, the total number 


of delay free feedback paths for complex wave digital filters is five, again 
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localized towards the source end. Had we designed the wave digital 
filter with no delay free path in port one this localization would 
have occurred towards the load end. Now with the above points in 
mind we can establish the following: 

(a) From the results of the large number of different filters 
studied in this chapter, it is apparent that wave digital filters do 
exhibit a high sensitivity to terminating resistances. For the case 
of filters with low terminating source resistances (i.e. R.-0 and 
ΒΞ) the rms error behavior of the wave digital filter is almost 
identical to direct digital filter derived from the analogue LC struc- 
ture, because of the localization of the delay free feedback paths 
towards the low resistance termination. For the case of the filters 
with high source resistance, the localization of the delay free feed- 
back path occurs towards the high resistance termination; thus the rms 
error, compared with that of the direct digital filter, is worse; this is due 
to the localization of the delay free feedback in the high sensitivity 
port of the filter. Над we designed the high resistive source 
impedance filters with no delay free path in port one, (thus directing 
the localization of the delay free feedback paths towards the load end) 
we would have obtained 2 better performance. 

(b) The rms error performance of simple wave digital filters were 
considerably better than complex wave digital filters for low source 
resistances (i.e. К 50, К 51). Thus rms error behavior can be explained 
as follows because there are fewer delay free feedback paths for the 
complex wave digital filters. This feedback occurs in fewer sections 
and the error performance gets worse. For high source resistance, i.e. 


R.-10, the rms error performance for all wave digital filters after 
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fifth or sixth bit becomes almost identical. This can be explained 
because of the high sensitivity of the wave digital filters to the 
terminating high source resistance. Thus the localized delay free feed- 
backs near the source increases the error performance for the simple 
wave digital filter while for the complex wave digital filter, with 
fewer feedbacks, the sensitivity to high source impedance drops.  There- 
fore the two algorithms tend to give nearly the same rms error per- 
formance. 

(c) In some cases especially for lower source impedances the 
ms error performance of reduced complex wave digital filter is slightly 
worse than that of complex wave digital filter. This is probably due 
to the fact that, when the truncation process is done on a large number 
of multipliers, the effect of truncation tends to be compensating, but 
if the number of multipliers are reduced, this equalization is less 
evident. 

(3) From the above discussion we can conclude that the ultimate 
ms error performance of wave digital filters tends to that of the 
directly digitized analogue resistively terminated LC filters. To achieve 
this ultimate performance in designing the wave digital filter, we have 
to use simple sections as much as possible, and also since terminating 
resistances are the most important we must direct the delay free port of 
Wave digital filters towards the low resistance termination. 

(4) If we wanted to make a band pass or band stop filter, which 
has LC tank circuits, the complex wave digital filter derived from 
these LC sections would not give the same ms error performance as that of the 
original conventional digitized LC filter. Moreover the algorithms LOr 


complex wave digital filters derived from LC tank circuits in the 
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literature [2], [3] are of the reduced parameter type complex wave 
digital filter, which furthermore have poorer error performance as 
can be seen from the graphs in some cases. 

In Chapter VII we discuss briefly how to design a simple wave 
digital filter algorithm for a given band pass or band stop specifi- 
cation. It is expected that these filters will have better error 
performance than the complex wave digital filters given in the 
literature [2] and [3]. 

(5) it is also interesting tonote that, after familiarization with 
wave digital filter theory, it is much easier to design a wave digital 
filter derived from analogue LC filters than to design a direct con- 
ventional digital filter both derived from the same analogue LC 
filter. Thus the wave digital filter is a practical design technique. 
Furthermore in order to reduce the rms error performance of wave 
digital filters we can even make use of a combination of complex wave 
digital and simple wave digital algorithms. For example when we have 
a high source impedance, we can make use of one or two complex sections 
near the high resistance end, followed by several simple sections later 


On. 
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VI. DERIVATION OF PARTIAL SENSITIVITY 
FUNCTIONS OF WAVE DIGITAL FILTERS 
A. INTRODUCTION 
The main intent of this chapter is to derive the sensitivity of 
wave digital filter with the aid of partial differentiation of the 
filter algorithm with respect to wave digital filter multiplier coeffi- 
cients, i.e. o's and ¢. Having found these sensitivity functions the 
Sensitivity of wave digital filter with respect to the original filter 
component values, i.e. L's, C's, Re, and Rr were found. 
Finally the behavior of these two sensitivity functions are investi- 
gated and compared with each other in the frequency domain and with 
the results obtained in Chapter V. 
B. DERIVATION OF THE SENSITIVITY FUNCTION OF WAVE DIGITAL 
Επι ας τὸ VARIATION IN MULTIPLIER COEFFICIENTS 
The general wave digital algorithm for the third order filter given 
in Fig. 6.2, with no delay free path in port two, is of the form of 


equation (6.1) 
a5, (2) bz, (2) 
f, (03,2) £2(07,2) (6.1) 


511 (2) $652 (2) 








where Ё,, t, and ες are 2x2 matrices of the general form o£ 
-1 _-1 
EI ο. 
E Ор l+y D 
£. (c, ,z) cl ee 
E. a5+832 της 
—— ley,z | 





Fig. 6.1. A general third order doubly terminated analogue low pass 
LC filter. 


Section 3 


(04) 





Fig. 6.2. The wave digital filter derived from the analogue LC filter of 
Fig. 6.1 with no delay free path on port two. 





where 05, 8;, ү; аге functions of σι. 
From Fig. 6.2 the transfer function of the filter in the time domain 
as depicted in Chapter IV and also explained in Appendix 1A with unity 


impulse input will be 


hat) = =>) 54, (2) = h(o,, 9,,---9) (6.2) 


Thus the transfer function of the filter in the frequency domain will be 


N 
Е > | -JwnT 
Ηίσι, с), ф) E b; (nT)e 
n=0 


and the sensitivity of the filter due to variations in 91 in the fre- 


quency domain is given by 


oH (c 


E») . 1$ > 
304 


32 ο ПТ 


Note that from equation (6.1), the matrices £5; f; are not functions of 


σι. Thus the derivative of equation (6.1) with respect to 91 will be 
0 90, (2) 
df, (94, z) 997 
= Ele Е, (0,2) f3 (03,2) (6.3) 
9011 (2) 3 90, (2) 
804 | |. TET. 304 


Note a5, (2) is input and is independent of σι also. 
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A closer look at equation (6.2) reveals that with unity impulse input, 


the time domain sensitivity function of the filter with respect to oi is 


ЭВ (91 ,0>,---$) _ 1% 9552 


901 2 904 


Thus to find the sensitivity function of the filter with respect to 91 


all we have to do is write the complete iterative equations of the 
filter from equation (6.1) and then differentiate them with respect to 


G Note that in doing so actually all the iterative equations derived 


1: 

from matrices Ғ,(0,,2), f,(0-,2) will not be effected and the only 

iterative equations being differentiated will be that of matrix 

£, (07,2). So in general we can find 22 , i.e. the sensitivity of the 
n 


wave digital filter with respect to all wave digital filter multiplier 


9Η 
дф 


coefficients ба" For the special case of from equation (6.2) we have 


902 


οπ(πτ) 1 
=> (bs, QT) + 9 s—— 


дф 





and the filter sensitivity function with respect to ф in the frequency 


domain will be 


N 
aH (o4 ,0,---9) ћ w a F (1+5) ‚ P52, -junT 
-. ἘΠΕ 

n=0 


The above procedure can best be illustrated with the aid of a simple 
example. 
Example 1 

Given the third order wave digital filter of Fig. 6.2 which is 


derived from the doubly terminated LC filter of Fig. 6.1; it is required 


EL 





to find the sensitivity function of the filter with respect to 9, in 
the frequency domain. The wave digital filter is with no delay free 


path in port two. 


Procedure 
The iterative equations of the wave digital filter of Fig. 6.2 in 
the time domain with the aid of Table 4.2b and with all initial condi- 


tions set to zero are: 


b,,(n) * a, (n)*a,4 (n-1)-a, , (n-1) *, (a, ; (n-1)-b,, (n-1)) 
E - 57) 

Boom) = 2,,(n-1)+6, (a, (m)*2,, (n-1)-a,,(n-1)-b,,(n-1)) 
pu) - 555() 

bz,(n) * a, (n)*a., (n-1)-a,, (n-1)*6; (a, (n-1)-b,, (n-1)) Би 
as (n) = b.,(n)- 9 

b; (n) * a4, (n)*os (a. (n)-2.,, (n) *a;, (n-1)-544 (n-1)) 

o 

ba, (n) = a,,(n)-a,, (n)+a,,(n-1)+0,(a,, (n)-b,, @-1) 

pre b) (n) 


971 (n) = a, (n) *o, (2, 5 (n) -a44 (n) *a,5 (n-1)-b44 (n-1)) 
Note that to find the transfer function of the filter in the time 
domain, the input 8,1 (n) will be 


1.0 =0 
= 
0 nz0 
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and the iterative equations after partial differentiation with respect 


to 91 are 
3b, > (п) 9a, , (n-1) да. „(п-1)  23b,.,(n-1) 
12 = = G Z Е Б JR 12 
304 997 ji aj (n 1) b, (n ОЕ p 304 ) 
Hou — 3550 (6.5) 
997 304 . 
>> (п) Е 9а, (п-1) ἽΝ αν . 32,4 (n-1) | да, (n-1) | de Ὁ 
ὅση 907 2 904 991 304 904 
etc. 


and in this way, we differentiate all the equations in equation (6.4) with 


mespect to 94 up to 


da, , (n) _ 3b», (n) 

304 да. 

ab oa (i Ja- (m 1) db,. (n-1) 
T E и ur _ |? 1 _ TI 

301 ат: (п) атт (п) +ау,0- 1 b, ms до, до, 


Thus the sensitivity function of the filter in the time domain from 


equation 6.2 is 


опт) _ 1, a 
904 У, 304 





and the sensitivity function of the filter in the frequency domain will 


be 








Q2 


N 
oH(o) _ dh(nT) _-jwnT 
да o D» σ £ 

1 В 1 


=0 


Jl. 














where T is the sampling interval, o is the frequency of input signal, 
and N is chosen large enough for transients due to the impulse response 


to decav. 


C. DERIVATION OF SENSITIVITY FUNCTION OF WAVE DIGITAL FILTER IN THE 
FREQUENCY DOMAIN WITH RESPECT TO ITS ORIGINAL COMPONENT VALUES, 
МЕ L's, C's, К, AND Rr 


As noted in Table 4.1 and 4.2 and mentioned in Appendix 1-A, there 
exists a one to one relationship between the wave digital multiplier 
coefficients, i.e. g's and the reflection coefficient $ of the wave 
digital filter, with the original filter component values, i.e. L's, 
es, Ro» and Rr. Thus in order to find the sensitivity of the wave 
digital filter algorithm with respect to original filter component 
values, the best way is to find the sensitivity of the wave digital 


filter with respect to wave digital filter multiplier coefficients, i.e. 
3H(01 ,0,---0) 


5С , etc. and from these, the sensitivity of the wave digital 


il 
filter with respect to original filter component values, i.e. 


SH(L, ,C5,L---- ,R., Rp) 
28 
illustrated with the aid of an example, the results of which will be 


, etc. can be found, This procedure can best be 


used later on to find the semi-logarithmic sensitivity (or percentage 
sensitivity; or normalized sensitivity) Of a seventh order 

low pass resistively terminated wave digital LC filter with respect to 
original components (i.e. L's, C's, Ro» and К). Later on these nor- 
malized sensitivity functions will be compared with that of the 
normalized sensitivity functions of the wave digital filter with respect 


to wave digital filter multiplier coefficients (i.e. o's and 9). 
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Example 2 


Given the sensitivity functions of the seventh order low pass wave 


digital filter of Fig. 6.4, designed with no delay free path in port two 


from the original LC structure of Fig. 6.3; it is required to find the 


sensitivity function of the wave digital filter with respect to original 


wave digital filter component values. 


Procedure 
Since the wave digital filter is of the type with no delay free 


path in port two, from Table 4.2 we have 


nen) (6.6) 


where Ry = Ro 


Also at this stage it is important to note that L and C compo- 


nents of the filter used in equation (6.6) and in the remainder of 


ΩΤ 
the example are predivided or premultiplied by the factor tan ce / 


w- İn order to take into account the effect of the sampling period 
and also the critical frequency of the digital filter, where tog is 
the critical frequency of the digital filter in rad/sec and ш 15 
the critical frequency of the analogue LC filter in rad/sec and T 
is the sampling time in sec. This factor is obtained by the multi- 


Q 
plication of the prewarping factor, i.e. 2 tan pos ως by the 
sampling factor of the νε components of the wave digital filter, 


| С 2S T 
Е. 1/2. The factor tan E reduces to tan E when the 


critical frequency of the original analogue filter is normalized at 


m 1 гаа/ѕес. 
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εις, 6.3. 


Bus. 6.4. 


Seventh order low pass analogue double resistively 
terminated LC filter. 





The wave digital filter derived from the seventh order 
analogue LC filter of Fig. 6.3 with no delay free path 
on port two. 





Thus from equation (6.6) 


сл 
| 

а 

сл 


2 Wak 
where 
| 
G. = — , 1 = 1,2,3 
Also 
2 ος е7 > C 
and 
R; = OR, 


Also from equation (6.7) and (6.8) we have 


с. С ; 
T 
G, = ——— = f(c.,C.) 
2 976, +6) j: 
thus 


2 2 
1 (c, G, +C.) (0, +R,C,) 


In the same way 


E se 
ЕЕЕ ТЕ ег. = АЫ ЫЗ, 
μ.ο — το 
9σ- R.L; 


(6.7) 


(6.8) 


(6.9) 


(6.10) 


(6.11) 


(6.12) 





σ ΝΜ. un ey) (6.13) 
4 σι σι στο”, $554 
Әс. GSC, (6.14) 
90. .— 
(04 6-*C,) 
E EE ua Г.) (6.15) 
5 R +L; o R. +L 5 
E RALc (ORTO) 
00, — 2 
4 (S¿Ry+Lo) 
σος οτε 7 m. 5 (6.17) 
6 Ὁ 55 6 ° 
9G ας 
a ο 56 (6.18) 
EL LS i 
pn mu RL” 500617) (6.19) 
90 | Rel- 
> (6.20) 
6 (o R ç+tL;) 
77 + R) 
Rp * Rg R + Go; 7 (6.21) 
a _ ΤΙ 
9σ 2 
7 (Rp *G797) 19:22 


Note that from equation (6.21) it can be seen that 


D Е, C5; Lz, ОВА Lo,R.,R, ) 
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Note also that the transfer function H of the filter is a function of 


ΓΞ (04, ΠΠ 7? $) (0725) 


and also 


КЕЕШ Соу т 


1? 25 39777) RR) (6,24) 


Thus if we have 2——, we can derive $ from (6.23) by writing all 
S 


partials of том етом (6.25), і.е. 
S 


MEME. E a о ως + —— ç — + 














IR, 304 OR, 30^ IR, 30 = IR, 
ООО О η) 3% он 9% 
A = ° Ри. ө + 
90 9 90 IR да oR 
4 S 5 S 6 S 
дс 
ЭН. 7 oH , 96 
307 "OR. " 36 R = 
S S 
aH aH aH aH 
= • ч—— + AL — + А, • <—— + А < 
|| 304 2 Др 3 да 4 99 4 
| aH νη oH . 9H 
И А. e EN b 30, + Ag ὅσα (6.26) 


to find A,, Ay, MT etc. we have from equation (6.6) with R, = R,- 


904 Ly (6.27) 
“л R, πα 


and from equation (6.25) and the chain rule 


Q2 


i да | да 
904 д 


Q2 
Q 


ΓΩ 
= 


> 
II 


| 


N 
Q2 
АЈ 
А 


ui 
(n 





Thus from equations (6.10) and (6.27) we have 





u 
A, = Е аа o (6.28) 
sp 
and in the same way 
90. na 
А, = s= = À, • 2 (6.29) 
- (97R7+L3) 
дс G-C 
4 Sea 
A, = Ξ Ἂς « ---- ΓΞ 6. 50 
4 R Š G.+C.)2 Ios 
т 
да RE 
А, = == А, Е 1 ; (6.31) 
+ 
: (σι τς) 
9σ Gel 
ker. — (6.32) 
ý (©_б6.+С.) 
5 5 0 
об RE 
А. = = = А, + sI Gass 
S (g R +L.) 
. |. L τ 
Ag = DR = An < ει v (6.34) 
S О, 


Thus from equations (6.26) to (5.34) we can find the sensitivity of the 

wave digital filter with respect to Ro = Ry To find the sensitivity 

of wave digital filter with respect to PE deer а from equation (6.24) 
1 


we can write 
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9H - ος , 991 , 3H 902 ан. μέν 904 
aL, 3c, 3L, 3e, db, до. 5η до, L) 
3H 205, ан , 206, ӘН , 307, 3Н (фф (655) 
3c. δη 80 3] 99; 3 9% 3L. 
aH aH 3H 3H 
Ее В В а l 
1 301 2 30^ до, 4 90 4 
+ В ° 9} | B ° oH + B : Er (6.36) 


Eo find Ву, Bo, Bz, etc. we have from equations (6.6), and (6.36) 
дс = 
B ^a - ies (6.37) 
1 (R, +L, ) 


and from equation (6.55) and the Chain Rule 


Thus from equations (6.10) and (6.37) 





2 
a. ΕΔ 
В, B, * E σσ (6.38) 
(9,6,+673% +) © С) 
and in the same way 
a L 
B, = T = B. `° c - (6.39) 
g 
l (o5R,+L-) 
9g GG 
By = = = SS (6.40) 
|| (o,G,+C,) 
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R,L 


5 4 5 
B τι en m 
5 aL 4 2 
Jt (o R +L.) 
B - . E. τ᾽ 
6 ο], 5 2 
1 (σςῶς”(ᾳ) 
п> Ж Reb 
B. UE ^ Ps ` E 
Dus 77 
B _ °$ = B Eoo Z 
8 | ο 2 
(o 7G *R) 


Thus from equations (6.3 


wave digital filter with 


To find the sensitivity of wave digital 





(6.41) 


(6.42) 


(6.43) 


(6.44) 


6) ta (6.44) we can find the sensitivity of the 


ВЕРЕ 6 L4. 














filter with respect to C 





2? 
i.e, S from equation (6.21) we can write 
2 
35 91, sq | 392 эн. 55, ан, 294 
902 304 9С, 90, 902 30. 902 90 y 9 2 
EEUU ен ε Ἡ 777 эн ә 
boxe эз ος 30. 3C. ὃσ ^ 30; (6.45) 
5 2 6 2 7 2 2 
Br 
aH aH aH aH aH aH 
О. +). 6 р < + 0 <—— 
aC, l 991 2 90, 5 90, 4 304 5 90, 
oH aH aH 
Е раар. (6.46) 
6 906 7 907 8 30. 
= 0. 


Note that the coefficient of the first term, L.e. Di 


From equation (6.9) 
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2 2 
Lo στ Y (6.47) 
2 9 2 (6202) 
and from equation (6.45) and the Chain Rule 
90. да, 30^ 
7s 7 36, " 3o, 30, δν 
Thus from equations (6.12) and (6.47) 
R.L; 
D; = D, : ра E 7 (6.49) 
a: 
and in the same way 
G.C, 
D, = D; r (6.6.6 97 (6.50) 
ENS 4 
RL. 
pp = (6.51) 
(o,R, +L.) 
Sat 
De = D- ° (бк)? (6.52) 
5 5 0 
ESL 
Па = De : AI, (6.53) 
(g R +L.) 
-2G 
E. con (6.54) 
8 7 Z 
(c-G*R) 


Thus from equations (6.46) to (6.54) we can find the sensitiyity of the 


wave digital filter with respect to С. 


To find the sensitivity of wave digital filter with respect to 


EC 


3) L. in the same way we have 


4» bs» Cg, Lo 
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and 


E EN 
(o,R,+L 


4 50, "53 











5) 


6566 


2 


RL, 


"267. 


᾽ (9565*С | 


СИЛЕР 2. 
(g R +L.) 


о 2 
(9-6. К) 


aH H aH 
Dt o Bor Eo* =— 
6 бе 7 99 8 90 
aH aH 
а. 
aH 
m 
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i955) 


(6.56) 


(6.57) 


(6.58) 


(6.59) 


(6.60) 


(6.61) 


(6702) 


(6.63) 


(6.64) 


(6.65) 


(6.66) 


and 


and 


and 


and 


σσ) - 
„Мз 
ов.) 
Gs, 


> (0.6, Br 


δε 


om 


К 
| -2G.R, 
ee ο 
С7*В 


-Rs 


E 2 
(Ro+L.) 


ους 


o Ge *Cc) 
nx 


ui Re ЕЕ 7) 


чү 
7 (g CR) 
ТТ. 
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(6 


(6 


(6. 


(6. 


(6. 


(6. 


ο 


(6. 


(Or 


(б 


(б. 


(O. 


(6. 


67) 


.68) 


.69) 


70) 


Lb 


p 


15) 


74) 


75) 


76) 


ΓΙ, 


78) 


79) 


80) 





ES 26.0, 
v = ——ə — (6.81) 


(c. G,*R )^ 

Thus we can find the sensitivity of wave digital filter with respect to 
the original filter component values. 
D. EXPERIMENTAL STUDY ON THE INTERNAL SENSITIVITY 

BEHAVIOR OF WAVE DIGITAL FILTERS 

In order to analyze the internal sensitivity behavior of wave digital 
filters with respect to the multiplier coefficients and compare them to 
the sensitivity of wave digital filters with respect to the original 
filter component values, in total nine different seventh order low pass 
filters were taken from the Handbook of Filter Synthesis [1] with source 
resistances varving from 0 to 10 ohms. The sensitivity behavior of these 
filters in the frequency domain with respect to both wave digital filter 
multiplier coefficients and the original filter component values were 


found using the procedures set in examples l and 2. These filters were 


m seventh order .5 db ripple low pass Chebyshev filter with R.=1.0 
ii) seventh order .1 db ripple low pass Chebyshev filter with R.-1.0 
lii) seventh order Butterworth low pass filter with R,=1.0 
iv) seventh order .5 db ripple low pass Chebyshev Ве ВОО 
v) seventh order .1 db ripple low pass Chebyshev filter with R =10.0 
vi) seventh order Butterworth low pass filter with R_=10.0 
vii) seventh order .5 db ripple low pass Chebyshev filter with R,=0.0 
viii) seventh order .1 db ripple low pass Chebyshev filter with R,=0.0 
ix) seventh order Butterworth low pass filter with τ 
Note that these filters are the same filters investigated in the sensi- 
tivity analysis of Chapter V. In order to make the comparison between 


the two sensitivity functions, i.e. sensitivity with respect to wave 
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digital filter multiplier coefficients and sensitivity with respect το 
wave digital filter original components more meaningful, the semi- 
logarithmic sensitivity function which is in fact normalized or percentage 


change in sensitivity is used. These normalized sensitivity functions 


are — : m and are plotted on the same graph in the frequency domain 
n 


MEETS. 6.5 to 6.11. 
Note that there are nine variables for the seventh order filter in 


terms of original components (i.e. Li; C»; Lz, р и. В. мата 


4 5? б” Т? SE 
R.) while there are eight variables for the filter with no delay free 


path in port two (i.e. 915 05,---0- and $), since the source reflection 


Coefficient, € is made equal to zero. Thus the curves of = - R_ are 
S 
plotted on a single graph. The remainder of eight variables are plotted 


on the same coordinates, i.e. S c L4 is plotted on the same coordinate 


aH aH . OH 
aS 377 σι and so on. Finally R R. 1s plotted with s= @ on the 


Same Eu: It is important to note that for the filters with 

Ro UU. tne normalized partial sensitivity <= - 91 cguals torero Se 
су + O when Ry = R. = 0. Also we note that РА this special case 

* lue "nveguation (6.26) all the coefficients of Sp are zero 


with Ry = Ro ^ 0. Thus for these reasons for when R, = 0 we have 
H Πα. 
plotted only 25 rather than =. ci and we also note that a is 
В i J 
Zero for all frequencies. 


E. ANALYSIS OF THE INTERNAL STRUCTURE OF THE WAVE DIGITAL FILTERS 
IN FREQUENCY DOMAIN USING FILTER PARTIAL SENSITIVITY FUNCTIONS 
Although all the filters under investigation are seventh order low 


pass, they were chosen to be as different from each other as possible, 
in type and source termination resistances. As mentioned earlier these 


filters are identical to those analyzed in Chapter V for sensitivity 
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in the frequency domain. Thus it is expected that the product of the 
sensitivity study of this chapter to conform to the bit 

truncation study of the Chapter V. It is important to emphasize that 
in this chapter only the internal sensitivity behavior of the simple 
wave digital filters is considered. 

Analysis of the normalized partial sensitivity function curves of 
Pres. 6-5, 15 reveal the following. 

1) Both sensitivity functions, i.e. sensitivity due to wave digital 
filter multiplier coefficients and sensitivity due to the original wave 
digital components peak at the critical frequency in all cases. This 
can be expected, since it is well known that the digital filters in 
general exhibit high sensitivities around the critical frequency of the 
filter. The critical frequency of the filters under investigation was 
1 radian/sec and this point has been marked by a vertical line on 
ate sensitivity graphs in Figs. 6.5, 13. 

2) In general the normalized sensitivity due to multiplier coeffi- 
cients are smaller than sensitivity due to original filter components 
for the first sections. Gradually this difference gets smaller and 


for some few cases in the last section it reverses. This is reasonable 


Since it can be expected that earlier multiplier coefficients will have 
lower sensitivity than the coefficients associated with the later sections. 
3) At low frequencies, the normalized sensitivity of the wave 
digital filters due to multiplier coefficients are slightly larger than 
the sensitivity of the filter to the original filter components. This 
effect reverses in the mid and high frequencies, 
4) The sensitivity function of the wave digital filter with respect 
to the source resistance increases with increasing source resistance, 


asvevedent from Figs. 4.5, 13. 
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5) The sensitivity with respect to terminating load resistance and 
reflection coefficient have exactly the same pattern as evident from 
equation (6.60). Apart from this fact it is important to note that both 
sensitivity functions tend to follow the frequency characteristic of 
the original filter. This can be seen from the graph 6.13 for the case 
ἃς - 10. As R. decreases the fluctuations in the sensitivity curve in- 
creases while following the same pattern. In fact for R9 these fluctua- 
tions increase to such an extent that at some frequencies the sensitivity 
becomes negligible. This phenomenon is also evident in the multiplier 
coefficients of the last few sections of the wave digital filter but 
pora lesser extent. 

6) In general the sensitivity of the wave digital filter due to 
variations in multiplier coefficients is lower than the sensitivity of 
the wave digital filter with respect to the original component values. 

7) Also it is worthwhile to note that the high sensitivities 
observed for high source resistances in wave digital filters algorithms 


with no delay free path in port two are in complete agreement with the 


results obtained in Chapter V. 
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VII. CONCLUSION 


A. INTRODUCTORY REMARKS 

It is important to note that in this thesis, whenever appropriate, 
detailed conclusions and discussions o£ the results are made in the 
chapters concerned. Thus it would be repetitious to state these con- 
clusions again. However for completeness, the highlights of important 
Conclusions based on experimental results are summarized here. It is 
also important to note that in the experimental studies of this thesis, 
in order to achieve reliable and accurate results, a large number of 
filters of different types with different termination source resistances 
were studied, and the conclusions made on the basis of collective 
Results. 
B. SUMMARY OF THE IMPORTANT RESULTS IN THE FREQUENCY 

DOMAIN BEHAVIOR OF THE WAVE DIGITAL FILTER 

1 - Digital filters, derived from doubly terminated LC analogue 
filters using the bilinear transform, have the lowest sensitiv- 
ity of frequency response for variation of the original L, C, 
Ro and R. parameters o£ the algorithms tested. 

2 - Wave digital filters, derived from doubly terminated LC analogue 
structures, if designed properly tend to achieve exactly the 
same low sensitivity as that of the above conventional digital 
filters. 

3 - Design of conventional digital filters from LC structures is 
relatively a tedious job, while designing the wave digital 


filters from LC structures is relatively simple. 
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4 - Wave digital filters exhibit high sensitivity to termination 
resistance values. Thus in the design of wave digital filters, 
in order to achieve a desired performance, the delay free loop 
should be chosen at the low impedance termination of the filter. 

5 - Due to the internal structure of the wave digital filter com- 
posed of sections with multiple LC resonant elements, these 
filters exhibit higher sensitivity than the same filter made 
up of sections with simple elements. 

6 - The sensitivities of the internal sections of the wave digital 
filter, both with respect to filter multiplier coefficients and 
original filter RLC components, tend to peak sharply at the 
critical frequency of the wave digital filter. 

7 - The sensitivities of the internal sections of the wave digital 
filter increases towards the load end. 

8 - For the seventh order low pass wave digital filters studied, 
the slope of the mms error due to quantization in the number of 
bits of the multiplier coefficients versus the number of bits 
ls approximately 3 db per bit. Also, interestingly, the slope 
of the rms error of the conventional seventh order digital 
filter is also approximately 3 db per bit. 

C. SUMMARY OF NEW THEORETICAL EXTENSIONS 

TO WAVE DIGITAL FILTER THEORY 

1 - Earlier researchers [1] and [2] have stated as a typical example 
that for the port two resistance of the wave digital section 
designed for series L with no delay free path in port two,one 
should use R, = R, + L. This formula does not take into account 


the effect of sampling time. In the theory developed in this 





thesis, following a derivation parallel to the one given in [1] 
and [2], the port two resistance for the said section is derived 
as К, = R. + E . This result is more general than the previous 
one, for which an inherent fixed sampling time of one second 
must always be assumed. 

2 - In the theory developed in this thesis for the design of wave 
digital subsections, only one algorithm is derived for a multiple 
LC resonant section. This result is applicable for both series 
or shunt elements. This new approach makes it possible to design 


several alternate wave digital filter algorithms, for a given 


moño ee LC filter structure, 


D. SUGGESTED FUTURE RESEARCH 

1 - In the analysis of this thesis we have established that the 
Simple wave digital filter, due to its inherent ladder-like 
Structure, exhibits lower sensitivities to multiplier truncation 
than the complex cascaded section wave digital filter. Thus the 
need for the design of simple wave digital filters, even for 
band pass or band stop applications, arises. In the literature 
all the present algorithms for band pass or band stop wave 
digital filters are of the complex cascaded multiple LC element 
type. However it is possible to design simple wave digital 
filter, even for band pass or band stop applications using 
simple sections, so as to reduce the sensitivity of the overall 
structure. The technique basically is as developed in Chapter IV. 
The suggested outline is as follows. The analogue band pass or 
band stop LC filter derived from the low pass analogue LC filter 


would have sections of the type shown in Fig. 7.1. The design 
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Fig. 7.1. Total combination of L and C in a two port network. 


o 


195 





of simple wave digital filter for Figs. 7.1a and 7.1d is no 
problem and they can be separated into two port simple sections 
as shown in Figs. 7.2a and 7.2d. Thus we have to discuss only 
the cases of Figs. 7.1b and 7.1c. To design a simple wave 
digital algorithm for Fig. 7.1b we must employ a three port 

wave flow network. A general three port wave flow network is 
shown in Fig. 7.3. Note that in this figure there is a feedback 
path from all inputs to all outputs. Fig.7.2b reveals that in 
order to design two simple wave digital sections we have to 
employ a three port signal flow graph. Either L or C of Fig. 7.2b 
can be adapted into a three port structure. Su EET Е 
as the three port element. It can be shown that in order to 
have a causal network, two of the three ports must have no delay- 
free feedback path. We can make port-two and port-three of 

Fig. 7.3 with no delay free feedback. In doing so we can find 
R, in terms of К, Rz and L. Also we can find Rz in terms of 
Ry» R, and L. Thus we have two equations and two unknowns and 
we can solve for R, and R, in terms of R; and L. Now we can 
cascade a two port network designed for the component C into 
port three, and also cascade the succeeding sections into port 
two. It must be emphasized that the two port network designed 
for element C has no termination resistance on its second port; 
thus it is open circuited and its reflection coefficient in port 
two will be equal to 1. To design a simple section wave 

digital filter for Fig. 7.1c, we have to use a somewhat differ- 
ent approach. With the concept of delay free feedback in mind 


and with reference to Fig. 7.4, we can design a two port network 
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a) two cascaded two ports b) three port element connected 
to the two port element 





C) two series two ports in d) two cascaded two ports 
shunt configuration 


Fig. 7.2. Separation of the complex two port networks shown in 
Fig. 7.1 into simple two or three port networks. 





Port 1 





Fig. 7.3. A three port signal flow graph, with all possible combinations 


of delay free inputs and outputs. Note that the delay free 
path in port two and port three are marked with *. 
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for the element L such that there is no delay free feedback 
from port two (input) to port one (output) for subelement L. 
From this constraint, we can find R5', with Ry! arbitrarily 
made equal to zero. The signal flow graph for the element C 
then 1s a normal two port signal flow graph with no delay free 
path in port one. It is fairly obvious that port two of the 
composite section will have a termination resistance of R = 
R, + Rs' where R,' is the port two resistance of the element C. 
In the analysis of the internal structure of Chapter VI, we 
found the sensitivities of the wave digital filter with respect 
to the original component values, i.e. L's, C's, Re, and R, . 

It is interesting, for a given filter, to differentiate the 
conventional digital filter transfer function with respect to 
components L, C, n and R. and compare these sensitivities 

to those of the wave digital filter obtained in Chapter VI. 

The sensitivity functions of the wave digital filter with re- 
spect to filter multiplier coefficients can be used to imple- 
ment an adaptive wave digital filter, by feeding back a weighted 
percentage of the sensitivities into the input as shown 
schematically in Fig. 7.5. The adaptivity can be done using 
any of the known optimization methods available in the 
literature, namely gradient optimization techniques, Fletcher- 


Powell optimization techniques, etc. 
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Fig. 7.5. A schematic diagram for the proposed adaptive wave digital 
rer. 
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APPENDIX 1 


A. DESIGN OF LOW PASS LC FILTER WITH THE GIVEN SPECIFICATION 
specification 
It is required to design a .5 db low pass Chebechev filter with 
load resistance Rr = 50 2 and source resistance К. - 100 Q, with 
critical frequency of 100 radian/sec. 
2. Пата 
From the Handbook of the Filter Synthesis by Zverev [1] for the 


given specification, the normalized values of L and C for R, » 1 Q, 


L 
making Ro = 2.0 Q as per Figure А,1 аге 


Be 2.0 ohms 

L, = ,4799 henries 
C, = ,3536 farads 
L. = 72.2726 nenries 
Cy = ./512 farads 
Le = 3.5532 henries 
Ce = ,9513 farads 
L- = 3.0640 henries 


3. Design of Wave Digital Filter 
Using this data the required wave digital filter with no delay 
free path in port two was designed using the table 4.2b. The schematic 
wave flow diagram of the wave digital filter is shown in Fig. A.2. 
Note that from Figure A.2 and equation (4.17) the unity input impulse 


Pesponse of the filter will be 
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Fig. A.l. Seventh order low pass analogue double resistively 
terminated LC filter. 





Fig. A.2. 


Seventh order low pass .5 db ripple Chebyschev wave digital 
filter designed with no delay free path on port two after the 


seventh order low pass filter of Fig. A.1. Note that only 
delay free signal paths are shown. 
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etl) = ha) = 555 b- (nT) (A.1) 


where T is the sampling period of the filter. We also note that nini is 
also the transfer function of the filter in the time domain. Thus the 


transfer function of the filter in the frequency domain will be 
Но) = > вет) е 2911 


Or 


ii 
М 


1+6 -} 
E (=) b,,(nT) +e Juni 
n=0 
N 
l* -junT 
= Le > b+, (nT) ʻe J (A. 2) 
n=0 


where 9 is the reflection coefficient of the filter and N was chosen large 
enough for transients to decay. Using these results the frequency trans- 
fer function and also the impulse response of required filter was pro- 
grammed in the computer. Sampling period of .01 sec was used, thus 
making the sampling frequency approximately six times the critical frequency, 
or three times the Nyquist frequency. 

The appropriate computer programs for filter transfer function, and 
also its impulse response are given at the end of this appendix together 
with the computer output results which are listed in Tables A.1 and A.2 


with their corresponding graphs in Figs. A.4 and A.S. 
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4. Design of the Conventional Digital Filter for Comparison Purposes 


The transfer function of the given filter of Fig. A.1 is of the 


form 


V 265) 
H(s) = v (A. 5) 


It can be shown that H(s) is of the form 


where from 


V, (s) 1 


H(s) "TS —F m (A.A) 
| Ьуз «bes^«b. s? «bs «b.s «b, s^*b, sb, 


the circuit analysis we find 


K CL, 
S 
(A.5) 
" CE C LCS ap * Ly g2) 
L 
= 
= ung сро ост (Е + Cy) 
L 
27 
+ L;C L C  (CoRs + RD 
oem "S m e Cu ue lan pic. 5, с) 
χο R 
ο ей 
- E O GERE + 1 = Dr 


- (C 





m" 
+ Сабо (СЕСЕ) + = C6) 


НО aL a; [iE Ж О Е р Е 
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NEN! 
b, = в sb) τ 


and 
— 
O Rr 


The equation (4.4) can be rewritten in a more familiar form of 





bs V,(s) Ку 
ΕΝ 5 (A.6) 
Y Е s +ass tass +a s +a s +a s +a s+a, 
AI 
where K, = — 
n b; 
ps 
6 Da 
/ 
b; 
ἃς суа 
/ 
etc. 


To find the filter scale factor we let s > 0. This leads to the value 


+ L 
175 WO Ci iS the filter's scale factor. 
L 


To find T(z) for the direct digital filter design, i.e. the 
digital filter transfer function, we can use equation (A.0) and the 


bilinear transform equation (2.4) to get 


T(z) = H(s) 
(A.7) 





or for simplicity we can factor H(s) into one first order and three second 


order sections as shown in equation (4.8) 
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H(s) = ΕΝ. ...:.-. . 


2 u 2 A.8 
(s +A, S+B,) (s +4,5+В,) (s *A-S*B.) (S*A,) i» 


Note that in effect we are going to have one first order filter cascaded 
With three second order cascaded filters. It is much easier to bilinear 
transform the subsections one at a time rather than bilinear transform 


the whole H(s). Thus each second order section becomes 








1 
52+А $+В 
1 1 -1 
2s 
s = = ( i) 
| 7 
ERE. ee 
Е Bees з _-2 
|: 1 
where 
н 1 
К. = A 
р LEE. 
τ Г Í 
4 
5 - 8 
T^ 
cic? 
| 4 2 
A 
и Е 
т2 Т jl 
υπ. 
Eo 2 
τό Τ 1 





1 
S d AA 
У 2 dE 
S = = 
D 











E МЕ 
m -1 
1+9 42 
where 
1 
ГЕ 
5. 2 
T^ 
а 
и 
4 2 
+ A, 
inus (A.7) and (A.8) lead to 
iE E И " 
ссе) 2 | a —— 0e пат 
Ë - - 2 -1 -2 -1 -2 EI Ww 
(1+0, 2 +312 ) (1+0,2 +8,2 ) (1+u,2 +852 ) 1+0 42 


where 


Note that T(z) is merely the transfer function of four cascaded first 
order and second order sections as shown in Figure À.5. Thus the appro- 


priate iterative equations are 


Vi(nT) - KV; (nuT) * 2KV; (nT-T) + KV. (nT-2T) - oV, (nT-T) - &,Vj (nT-2T) 
Муса) = V (m) + T-T) + V,(nT-21) - a V,(nT-T) - 8,V,(nT-2T) 
Мат) = Ү,(аТ) + V T-T) + V5(nT-2T) - o,V,(nT-T) - 84V. (nT-2T) 
ШО) = ντα) — Vg(nT-T) - a,Vg (nT-T) (A.9) 


where T is the sampling time of the filter and with initial conditions 
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T4(z) 
2nd order 2nd order 2nd order lst order 
section section section section 


Fig. A.3. Seventh order digitized and cascaded filter corresponding to 
seventh order analogue filter of Fig. А.1. 
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set to zero, i.e. 


Tee) = Ü = 5 
V (I) =0 
ро 


Also we note that for impulse response 
1 for n=0 
Vig) * Vip T) - 


0 for n#0 


Thus the unity input impulse response of the filter from equation (A.9) 
will be У (пт) and we note that this is also the filter transfer function 


in the time domain, i.e. 
m mI) = h(nT) 


Also the filter transfer function in the frequency domain can be found 
from equation (A.2). Using these results the transfer function and also 
unity impulse response of filter was programmed in the computer and 
computed using the same sampling frequency as that of the wave digital 
Fer of Section 3. 

The simulation computer programs are given at the end of this 
appendix and their corresponding outputs are given in Tables A.l, A.2 
together with the results obtained for the same filter using wave digital 
filter design for comparison purposes. Note that no graphs are given 
for this simulation, since they were exactly the same as the wave digital 


filter graphs. 
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Fig. A.4. Sraph of the transfer function of the wave digital filter 
with the given specification in Section 1. 
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Fig. A.5. Graph of the unity impulse response of wave digital filter 
with the given specification in Section l. 
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Graphx, subroutine is given in Computer Program No. 1. 
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* 
E. Computer Program No. 4. Program for transfer function of the con- 
ventional cascaded digital filter with the given specification, 


ς **X*3 FREQUENCY RESPONCE жж жж 
ς CONNVENSIONAL DIGITAL FILTER DESIGN 
Е 
[ОЕ στ REAL?3 (A-H;C-Z) 
CO δι, EX*l16 wT;,Z | Y YH 
DIVERSICN ОАТА(210, 3) 
, DIMENSION (а). уст 
C 
C CCMPONENT VALUES 
Ç »35335***.503 CHEBECHEV LOW-PASS FILTER WITH RS-2.0 sso dos 
€ INDUCTANCE ANC CAPACITANCE VALUES IN FENRIES,ANL FARADS 
€ NORMALIZED TO CRITICAL FREQUENCY OF 1 RAD/SEC АТ 3 08 ECINT 
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О 
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Ді 524. 488506 
Co=.8137D0 
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18 РОМА" (10Х ,' ΤΕΙ, orp OAs. 'a S=", £12.53 //) 
С πη GF ide COEFFICIENTS UF 
C cM Ce eU eU ARCA Oecum caede eee 
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Е ο ОО О К лс Оа РЕМ PATA 15 COEF 
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ar о 2 
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TS MAL TPA S4AL SAL 7) /2L+2S%(C2+C 4+C €) 


κ 
Graphx, subroutine is given in Computer Program No. 1. 
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Program to calculate the rms error due to 
truncation in the number of bits of wave 
digital filter multiplier coefficients. 


APPENDIX 2 
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A - Computer Program No. 5. 
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Subroutine Graphx is given in Computer Program No. 1 of Appendix 1. 
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C - Computer Program No. 


Program to calculate the rms error due to 
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truncation in the number of bits of reduced 
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D - Computer Program No. 8. Program to calculate the rms error due to 
truncation in the number of bits of the 
RLC components of the conventional direct 
digital filter. 
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“Subroutine Graphx is given in Computer Program No. 1 of Appendix 1 
and Function Trunc is given in Computer Program No. 5. 
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CALL GRAPHX(DAT A3,985s 4HFR EQ; 4HMAG. ) 
WRISREN S5 4.) 
HEERES) _ 

P OR AT(23X, 'FREQ',11X,'(D9(0/P)/D(SIGMA4)]*SIGMA4',4X,'[D(0/7P)/D(L4 
WR ITE(6,20) ((DATAG(N+M) ¿M=1,3) ,N=1,9€) 

CALL GRAPHX (DAT 24,98  #HFREQ, #hMAG. ) 
WR ITE(6,54) 
ТЕС БЭ) 

S ou eR 'FREC' , 11X,'(0(0/P)/0( SIGMA5))*SIGMAS',4X, ' ( DCC/ P)/0(L5 

BUE o NIRE AS D. Mzl,31)5,N-1,98) 
CALL GRAPHX(DATAS,98, GáHFR EQ ,6HMÁG.) 
WRITE(6,54) 

AR ITE(6,60) 

«ΓΑΙΑ (20X, 'FREQ',11X '(D(0/P)/D(SIGMA6))*SIGMA€',4X,'(D(0/P)/D(L6 

42 ITE(6720) CiDATAS(N y) 4 M= 1, 3) NS 1, €) 

CALL GRÁPHX (DAT A6,93, 4HF REO, +HMAG. 

WE ITE(6,54) 

BS Ал) f 

,FORMATL20X, TERE Gt 11 X,  10(0/0) /DU SIGNAT) )*S IGVAT! Ao! DDEC/P 1/0007 
WRITE(6,20) ((UDATAT (NM) 5V21,3) ,N21,98) 

CALL GRAPHX(UDATAT,98, 4HFREQ,AHMAG. ) 
HEINE 554.) 
WR ITE(6,62) 

,EORTAT (23. ! FREQ! 4L1X o! (0(0/P 1/0( PHT} ) #ABS( PHT) !, 4X, ! LOCO /P) /DIRL) 
WRITELO6:Z0) ((DAT88IN,V) ,V2 1, 3) ,N21,9€) | 
GAME APHX (DAT 46,398 » 4HERES, SHMAG.) 

FORM (20% , El2.5 +10X,E12.5 ,10X,E12.5) 
s TOP 
ENO 
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